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Abstract 

The non-equilibrium critical-layer analysis of a system of frequency-detuned resonant- 
triads is presented- In this part of the analysis, the system of partial differential critical- 
layer equations derived in Part 1 is solved analytically to yield the amplitude equations 
which are analyzed using a combination of asymptotic and numerical methods. Numer- 
ical solutions of the inviscid non-equilibrium oblique-mode amplitude equations show 
that the frequency-detuned self-interaction enhances the growth of the lower-frequency 
oblique modes more than the higher-frequency ones. All amplitudes become singular 
at the same finite downstream position. The frequency detuning delays the occurrence 
of the singularity. The spanwise-periodic mean-flow distortion and low-frequency non- 
linear modes are generated by the critical-layer interaction between frequency-detuned 
oblique modes. The nonlinear mean flow and higher harmonics as well as the primary 
instabilities become as large as the base mean flow in the inviscid wall layer in the down- 
stream region where the distance from the singularity is of the order of the wavelength 
scale. 
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1, Introduction and summary of Part 1 


The non-equilibrium critical-layer analysis of a system of frequency- detuned resonant- 
triads in boundary layers with and without mean pressure gradient was obtained in Part 
1 (Lee 1998b) 1 . Each resonant-triad is composed of a fundamental plane wave and a 
pair of subharmonic oblique modes. It is shown that the resonant-triads can nonlinearly 
interact within the common critical layer when their frequencies are different by a factor 
whose magnitude is of the order of the growth rates multiplied by the wavenumbers. 

If the (fundamental) frequencies of the resonant-triads are sufficiently detuned, the 
plane and oblique modes of each resonant-triad first nonlinearly interact between themselves. 
Their growth is not affected by the existence of the other resonant-triads and the growth 
rates are increased due to the parametric-resonance and self-interaction effects (of a single 
resonant-triad interaction as shown by Goldstein & Lee 1992). The frequency-detuned 
resonant-triads start to interact between themselves at the downstream position where the 
instability growth rates become large enough to be equal to the magnitude of the phase 
speed differences (of the plane waves). Eventually, as the magnitude of the growth rates 
approaches that of the wavenumbers, most of the unstable waves whose scaled Strouhal 
numbers are different by nearly 0(1) can nonlinearly interact. In this downstream region, 
the frequency range in which the resonant-triads can nonlinearly interact expands to cover 
the entire range of linearly unstable waves. 

The results of experimental study by Corke & Gruber (1996) (an extensive survey on 


3 f appeared in parentheses in the last term in (3.1) of Part 1 must be replaced by iflcjYcj in front of 
the square brackets in (5.6) of Part 1 must be replaced by i Y CJ , and in (7.60) of Part 1 must be replaced 
by 4 0) . 
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resonant-triad interaction can be found in Kachanov 1994) and direct numerical simulations 
by Liu k Maslowe (1998) in Falkner-Skan boundary layers are in good agreement with the 
resonant-triad theory of Goldstein & Lee (1992). The parametric- resonance-, self- and 
mutual-interactions and back-reaction are identified and shown to play major role in the 
transition process, especially in the adverse-pressure-gradient boundary layer. 

The system of partial differential critical-layer equations can be solved both analytically 
(Goldstein k Choi 1989; Goldstein k Lee 1992, 1993; Wu 1992, 1995; Wu, Lee k Cowley 
1993) and numerically (Lee 1997a). In this part of the study, the system of critical-layer 
equations given in §7 of Part 1 will be solved analytically and the amplitude equations 
(without the back-reaction term in the plane-wave amplitude equation) will be presented 
along with numerical solutions. The analytical solutions of the critical-layer equations are 
obtained in §2. The frequency-detuned non-equilibrium amplitude equations are presented 
in §3 and their large frequency- detuning limit and asymptotic singular solutions are given 
in §4 and §5, respectively. The viscous limit of the amplitude equations is obtained in 
§6. The velocity jumps across the critical layer for the nonlinearly-generated low-frequency 
modes are given in §7. Numerical solutions of the frequency-detuned quasi-equilibrium 
(i.e. viscous-limit) amplitude equations are presented in §8. The frequency-detuned non- 
equilibrium amplitude equations of the oblique modes are numerically solved in §9. The 
concluding remarks are given in §10. 

In the rest of this section we will give a brief summary of Part 1 to allow the reader to 
follow along without consulting Part 1 too often. 

In the non-equilibrium critical layer the mean convection effect balances with the 
growth and viscous effects. If we introduce the wavenumber parameter a and the local- 
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growth-rate parameter a r to characterize the small wavenumber and the ratio of the small 
local growth rate to the wavenumber, respectively, the generalized scaling of Lee (1997a) 
can be written as (see also (1.2.2) - (1.2.4) where henceforth an I at the beginning of an 
equation number will indicate an equation given in Part 1), for — J ^ j ^ /, 

ol 3 = a[a 3 + 0{a r )\, Cj = a[c 3 + 0(a r )}, c 3 = <r[c 3 +0(a r )] 7 (3 3 = afi 3 , y CJ = aY CJ , (1.1) 

Xl = C7 r + 1 X, (1*2) 

H — a T ~ l jl for 1 ^ r ^ 3 and fi = 0(1) for 0 < r ^ 1 , (1.3) 

where c 3 and c 3 are the (nearly equal) phase velocities of the two-dimensional and oblique 
modes respectively, a 3 is the streamwise wavenumber of the plane wave, j3 3 represents the 
spanwise wavenumber of the oblique mode, y CJ is the critical level where the base mean-flow 
velocity is equal to the real part of the phase velocity, a j7 c J7 fi 3 and Y CJ are order-one real 
constants and fi is the normalized mean pressure gradient. The subscript j is used to denote 
the quantities of the jth resonant triad. However, the subscript 0 for the quantities of the 
‘reference’ Oth resonant-triad will be omitted for notational simplicity. 

The local-growth-rate exponent r is a function of the streamwise coordinate since the 
instability growth rates are varying. The frequency-detuned non-equilibrium analysis of 
this paper is valid for any value of r. Each analysis for a specific value of r covers different 
stage of the non-equilibrium critical-layer evolution (see table 1 of Lee 1997a). A composite 
solution could be obtained to cover a series of nonlinear stages. 

The jth resonant triad is composed of a single two-dimensional mode and a pair of 
subharmonic oblique modes. The phase speed of the plane wave is assumed to be nearly 
equal to that of the oblique modes of the same resonant-triad. This resonance condition is 


: 1 I 
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satisfied when the propagation angle of the oblique mode 9 defined as (also as (1.5.29)) 

9 = cos -1 (1-4) 

is about 7r/3. 

The nondimensional frequency or unsealed Strouhal number s 3 of the fundamental two- 
dimensional mode of the jth resonant-triad is equal to a 3 c 3 and the scaled Strouhal number 
s 3 is given by (also by (1.2.5)) 

s 3 — a 2 s 3 = a 2 a 3 c y (1.5) 

The difference between the scaled Strouhal numbers of the jth and the ‘reference’ Oth 
resonant-triads is given by (also by (1.2.6)) (Lee 1997b, 1998a) 

Sj = 5(1 + cr r jx) or s 3 = 5(1 + jx) for (1*6) 

where 

X = X' C 1 - 7 ) 

It was shown by (1.5.22) that the scaled phase velocity of the jth resonant-triad, c 3 is different 
from that of the Oth resonant- triad, c by 0(cr r ), which is the additional resonance condition 
that is required for the nonlinear interaction between the frequency-detuned resonant-triads 
to occur in the common critical layer. 

In the non- equilibrium critical-layer analysis, the Reynolds number is scaled as (also 
as (1.2.8)) 

A = l/(a 3r+4 R A ), (1-8) 

and the amplitudes are scaled as (also as (1.2.9)), for all j, 

€U = tf 4r+ \ hd = *02 = tf 3r+1 , (1.9) 
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where €2 d 7 $3 d and 602 are the amplitude scalings of the plane wave, oblique modes and 
nonlinearly-generated low-frequency (spanwise-periodic) modes. 

The scalings given above are for the long-wavelength small-growth-rate instabilities 
in boundary layers, for example, in the high-Reynolds-number adverse-pressure-gradient 
boundary layer or in the downstream non- equilibrium stage of the Blasius boundary layer. 
Their frequencies are of the order of the wavenumber scaling squared and their critical layers 
are distinct from the viscous wall layer. 

The multi-layer structure for the non-equilibrium critical-layer analysis was given in 
figure 1 of Part 1. The unsteady flows in the main boundary layer ( y = 0(1)), the inviscid 
wall layer ( y = O(cr)), and the viscous Stokes layer 2 (y = 0(a 3r / 2 + 1 )) are governed by linear 
dynamics. The viscous Stokes-layer effect was included in order to make the analysis valid 
in the viscous limit. The solutions in the inviscid wall layer become singular at the critical 
level, therefore, they have to be rescaled in the critical layer. The nonlinear interaction 
between instability waves of the frequency- detuned resonant-triads first occurs within the 
common critical layer whose thickness is of 0(a r+1 ). 

The streamwise velocity in the main boundary layer is given by (also by (1.3.1)) (Lee 
1997b, 1998a) 

J j 

u = U + e 2 d 52 R e Bji x i)$jy(y’ x i) elXj + hd 52 cos Z 3 

3=-J J=~J 

+ <$02 52 £ ^ e Uo,2\j, t(y, £l> Z, ti) + . . . , (1-10) 

J l 

where 

Xj = (TOtjX — a 2 Sjt, Zj = afijZ, ( 1 . 11 ) 

2 y = 0(<r 3r/2+1 A 1/2 ) when A < 1. 


M II 
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U is the base mean flow velocity, A } and B 2 denote the oblique and plane wave amplitudes, 
respectively, and t\ is defined by (1.5.24). 

The jump equations are obtained from the requirement that the velocity jumps across 
the critical layer which are calculated from the critical-layer solutions are equal to those 
calculated from the external solutions. They must be solved with the system of critical-layer 
equations to determine the unknown instability amplitudes. 

The system of partial differential critical-layer equations along with the transverse 
boundary conditions and the jump equations are presented in §7 of Part 1. They are 
normalized in such a way that their nonlinear growth parts are free from any mean-flow- 
dependent parameter apart from A. The normalized variables in the critical-layer equations 
that are relevant to this part of the paper are, as in (1.3.19), (1.5.26) and (1.7.1) - (1.7.4), 


— K 2 0 ) i A 


— 3 1 X — - X) a lM 

CiC 6 KT W 


(2) _ vY c a r 3 _ r A y \ 2 

1M-£ T 3 ? ^ ° \2 o c ) ' 


( 1 . 12 ) 


A,/V lv (s) = (yc ^j B 2 /V 2;2 j (x) = Eehx^/2^ (1.13) 


where M = Sir YJ3 2 / (r^ac 3 ) as defined by (1.7.5), V au (x) = exp |i a(X 0 + r) 0 x ) - ijxx as 


given by (1.7.6), k is a normalization parameter which can be chosen arbitrarily, x 0 , X Q 
and T) 0 are the coordinate origin shifts (see (1.7.58) and (1.7.62)), t 0 is the scaled Blasius 
skin friction, t w denotes the total wall-shear stress (that is the sum of the Blasius skin 
friction and the correction due to the mean pressure gradient) as in (1.3.15) and (1.3.16) 


and x = X + 0(a) as in (1.5.23). 
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2. Solutions of the critical- layer equations 

The analytical solutions of the partial differential critical-layer equations given in §7 of 
Part 1 will be obtained in this section. 

As in the previous studies (Wu et al. 1993; Goldstein & Lee 1993; Wu 1995; Lee 1997a), 
the Iburier transform with respect to 77 will be used, 

/ co 1 roo 

e~‘ kr, Q(x, T])dT), Q(x,r 7 )= — / e lkv Q{x,k)dk, (2.1) 

- co J —co 

where Q = ^{Q} is the Fourier transform of Q. The velocity jump across the critical layer 
can be obtained from the Fourier transform by putting k = 0, 

/ oo 

Q(x,r})dr]. (2.2) 

-OC 

In order to simplify the presentation we will put 

£1 = it — (2 = x i - h s = t sin 2 0, (2.3) 

and 

t*(Ci;Ca) = exp [±A (|mCi + C 2 ) C?] • (2.4) 

2.1. The leading-order equation 

The solution of (1.7.8) which satisfies the transverse boundary condition (1.7.9) can be 
obtained by taking the Fourier transform of (1.7. 8 ), solving the resulting equation and then 
taking the inverse Fourier transform of <5^(ar,A;). The solution is 

Qj 1 ) = e~ l£T, I 0j (x, 7 ?), (2.5) 

where 

I n} (x,T})= j izie^^ffCijOXJSA^zi). (2.6) 

J — oo 
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2.2. The second-order equations 


In order to simplify the presentation we will only show the Fourier transforms of the solutions 
of (1.7.10), (1.7.11) and (1.7.14) which satisfy the boundary conditions (1.7.20) and (1.7.21). 
They are 

?{w$i i*E r dx 1 J J j(x,x 1 ,k), (2.7) 

?{w { 0 i, e , U$A (x, k) = 7T f dx x {i, — (2Ci + k )} [J jA {x, x l7 k ) + *i, *)] , 

( 2 . 8 ) 

r U^l h(r) } ( x , fc) = i7re Afc3/6 f dk{- i, fc} [e Afc3/6 fctf(-fc)Aj_*(x a )A,(x 0 + *) 

-2 h s e~ yk3 / e r dk i r dk 2 G-{k)A 3 . t {x a + ki)Ae(x a + * 2 )| , (2.9) 

J — oo J —oo 


where 


J^(x,Xi,fc) = Sl 1 (k;-Ci)kH(k)A J+ t(x 1 )Af(x l - k ), 

/ OO r OO 

■OO v/ — oo 

a*(Jfe) s e^ +fc i)/ 3 ^(-fex)^(-)fc 2 )^(ifc - Jfcx ± * 2 ), 


i a = x + (k - k)/ 2, 


( 2 . 10 ) 

( 2 . 11 ) 

( 2 . 12 ) 

(2.13) 


the summation notation is defined in (1.6.20) and the asterisk denotes the complex conju- 
gate. The step function H(k) is defined by (Butkov 1968) 


H{k ) = 1 if k > 0; H(k ) = 1/2 if k — 0; H(k) = 0 if k < 0, 


(2.14) 


and 6(k) is the delta function (Lighthill 1960; Butkov 1968). 
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2.3. The third-order equations for the oblique-mode velocity jump 


In order to obtain the velocity jump across the critical layer for the oblique modes, we need 
to solve the equations (1.7.22) - (1.7.30) (but only the second part which involves ^ >TO ]) 

with the boundary conditions (1.7.34). 

The linear growth term in the oblique-mode amplitude equation can be obtained by 
solving (1.7.22), 

? { y i,iyU} ’ k ) = 4in a[ 2 M eXk3/3H (- k ) A J (* + k ) » ( 2 - 15 ) 


and the linear velocity jump becomes, by putting k = 0 (see (2.2)), 


wM*)- (2-16) 


Using (2.5) and (1.7.13) we can show that the solution of (1.7.23) becomes 

^{nS^}(^^) = -Mcos 2 ^) /; dke l ^~ 2p ^ 3 k 2 H(k)B e (x + k- k)A}_ 3 (i + k- 2k), 

(2.17) 

and the velocity jump for the parametric-resonance term is 

= °) = ~ 47r ( cos2 0) / dari4“(Ci;0)Ci-Br(xi)A^_ J (2a: 1 - x). (2.18) 

^ ' J —oo 

The nonlinear part of the oblique-mode velocity jump is obtained by solving (1.7.24) 
to (1.7.30). It is easy to show from (1.7.24) and (1.7.29) that 

T KiL} (*, k = 0) = T (*, k = 0) = 0. (2.19) 

The solution of (1.7.25) can be obtained by using (2.7) and (1.7.10), 

?{vi\ ?L_}(x,fc = 0) = ^ f d Xl r dx 2 K^A^(x 1 )A e+m (x 2 )A* m (x 1 + x 2 - x), 

m: — l 4—00 J -<*> 

( 2 . 20 ) 
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where 


tf (3c) = £f(Ci;C 2 )Ci- ( 2 - 21 ) 

From (1.7.26) along with (2.5), (2.8) and (1.7.11), we can show that 

T {v/32^} (z, k = 0) = 7T dx a jT ds 2 Af 0 A J . e (x 1 )A e+ m(x 2 )Al i {x 1 +x 2 -x) 

+V& d \ ( 2 . 22 ) 


where 


(2.23) 


hf d) = 4“(Ci;C 2 )Ci [(Ci + C 2 X 1 + 2h s j^ 2 dC3^"(C 3 ; CiXCx + C 2 - Cs) 

/ x rx\ /‘oo roc 

dxi dx 2 dr) dk 3 g^ [J7* >m (xi, x 2 , * 3 ) + £r,m(xi, * 2 , £ 3 )] , (2.24) 

-OO J — OO «/ — OO J — OG 

g± = h s e -«* ±k 3 >”£f (Cl? 0)j£ [(*3 + 2C 2 )/o ( ^)(xi, 17) + /^(xi, 77 )] , (2.25) 

with (2.6), (2.10) and (2.11). Integrating by parts and changing the order of integrations 


using 


/ X tX\ rx tx 

dx 1 / dx 2 = / dx 2 / <fci, 

-OO J — OC •/— OO */X2 


(2.26) 


(2.24) can be rewritten as 


y( 3c 0 = —2rrh s f dxi j dx 2 K\] d J A t+m (xi)Aj- ( (x 2 )Al l (xi + x 2 - x) 

J — OG J — OO 

A J _Xx 1 )A, +m (x 2 )A^(a: 1 + x 2 - x)] , (2.27) 


where 


4/? = a c (Ci,0,0; 1) + 0}((i,O,O; 1), A'g? = a c (Ci,C2,0; 1)+ #( 6 X 2 , 0; 1), (2.28) 
£ c (Ca, Cfr, Co a) = dCs^e {a| 0 [Ca + 2C 2 + 3(C 3 - Cc)] + (l + 2A| 0 $) C 2 Cc} , (2.29) 
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, rC 1 

(Ca?C6>Cc|^0 = 1 d^g e 

Jo 

l Q d(4 {a£f (ic, id)ia + s 2 (£ 4 ] id)ic 



[l + 2A4(4 - Ca - Cfc) 2 ] } , 

(2.30) 

_ ^ — ^ |^|Cl+C2^<?4*^3C2+Cl^ (C2-C6) 2 + (Ci+C2“C6)(C3+2<2-2C6)C3+ ^Ci+2C2+ 3 C3^) C?J 

(2.31) 

0 = Ci ~ C 3 , 6 = Cl + <2 + C3, 

4 = Cfe + Cs ~ C 4 ? 4 = Ca + C 2 - Cb- 

(2.32) 


The solutions of (1.7.27) and (1.7.28) can be given as, using (2.5), (2.9) and (1.7.14), 
T { v i*'tl™nr,} (x,k = 0) = IT f dxj I dx 2 f dx 3 e- x ^ /3 CiA m (x 3 ) 

K J ^ — CO J — CO J —co 

C3^/(0; l)At~ m (x2)A}_ J (x2 + X3 — x) + 2 h s f d2 4 e~^ 3+ ^ c V 3 £7T (£ e ; — 1) 

«/ — oo 

A^ m (x4)AJ_j(a;3 + x 4 - x)j , (2.33) 

where ( 3 = x 2 - rr 3) £ 4 = x 3 - x 4 , the minus superscript in Qj is for the (3e) component 
and the plus one is for the (3/) component and we have put 

Gf( Ca; a) = 2 h s e~ l & 3 f+t> 2 ® /6 {if Li; 2 ± l) , (2.34) 

le = Cs + C4, i/S Cs + Ca, 4=0 + 2C2, 4 = 4 + Cl- (2.35) 

Using (2.26) and integrating by parts, (2.33) can be rewritten as 

?{ V Mrnr,v} (*> * = °) = 2 « h * f d '** H 

J —OO J — oo 1 

A*t-,{x i + x 2 - x) + /sr 6 (3e ' 3/) ^_ m (xi)A m (x 2 )^_ J (x 1 + x 2 - x)] , (2.36) 

where 

Al 3e ’ 3/) = 2 hs f Q l d( 3 J^ 3 dQg g [4+ (Co; C 2 XC 3 - 0) (1 + 2A44 2 ) + 0 + (Cs - C 4 ;C 2 )4] , 

(2.37) 

: 1 I! 
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(2.38) 


K (3e,3f) = K (3e,3f) + £ ^ (l + 2 Af tt g) T I.Ci] , 

<? 3 = e -2AKa+C 2 )C 2 C3^-( C i; C2 )^-( C2 - Ci)^-(C 3 ; 2Ci + 3C 2 ), (2.39) 

the (-) and (+) signs of (qp) in (2.37) and (2.38) are for the (3e) and (3f) components, 
respectively, and £ a and & are defined in (2.32). 

We can write the velocity jump obtained from (1.7.30), along with (2.5), (2.8) and 
(1.7.11), as 

fWfvilrnnv} = °) = 2 * h ° [ X dx 1 r dx 2 Kf h) A 3+( (x 1 )A m (x 2 )A* e+m (x 1 + x 2 - x) 

+yW, (2.40) 

where 

K? h) = 4“(Ci; C 2 ) [ G d( 3 i;(( 3 - Ci)Ci(C 2 - Cs), (2.41) 

Jo 

/ r rri roc roo 

dx a / d* 2 / dr) / ( 2 - 4 2) 

-OO •/ — OO j —OO j — OO 

and and <7+ are defined in (2.11) and (2.25). The above equation (2.42) can be rewritten 
as, using (2.26) and integrating by parts, 

V (3/l) _ _2 t T h s j dx x f dx 2 \Kf£ Am(xi)A 3+ t(x 2 )A* (+m (xx +x 2 -x) 

J — 00 J — 00 L 

+K^ ) A 3+e (x 1 )A m (x 2 )A* t+m (x 1 +x 2 - *)] , (2.43) 

where 

K { na = #(Ci,0,0;-l), Kf$ = aj-(Ci,C2,0;-l), (2.44) 

and Of is defined in (2.30). 
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2.4. Other third-order equations 


In order to obtain the mutual-interaction term in the plane-wave amplitude equation, we 
need to know the following solutions of the other third-order equations. 

The solution of the first part of (1.7.23), along with (2.5) and (1.7.13), which satisfies 
the boundary condition (1.7.34) becomes 


T (z, k ) = — 27re Xk3 / 3 + k)H(-k) — i(cos 2 0) j dke 2X ^ 3 ^ 3 kH(k) 

Bt(x + k — k)A*£_ 3 {x + k — 2fc)j . (2.45) 


Equation (1.7.37) subject to (1.7.43) is solved with (2.5) and (1.7.13), 


T 


{MS,/. y 3 ( SS J (*»*)“ I^Ccos 2 6) f dke l{k3+2k3)/9 {i, 2k) kH{-k ) 

B e (x b )A } -e(xb + k ), (2.46) 


where 

x b = x + (k - k)/ 3, (2.47) 

and we have used the fact / = 0 which is obtained by matching the solution of (1.7.42) 
with the outer solution. 


2.5. The fourth-order equations for the plane-wave velocity jump 

The velocity jump for the plane wave can be obtained by solving (1.7.44) to (1.7. 51) with 
the transverse boundary condition (1.7.54). The velocity jumps which produce the linear 
and mutual-interaction terms in the plane-wave amplitude equation will be presented. The 
back-reaction term in the plane- wave amplitude equation that is quartic in the oblique-mode 
amplitudes (Wu 1995) will be considered in future. 


The linear velocity jump is obtained from (1.7.44), 

?{U™} (*, * = 0) [= ujgkdn] = 2i xe^Bjix). (2.48) 

The mutual-interaction term is determined by considering the components, anc ^ 

U 2 fil,e n in (1-7-45) and in (1.7.49). The velocity jump by the U^ jtTI component 

becomes, using (2.7), (1.7.10) and (1.7.13), 

^K ( tL}(*’* = 0) = 87r £ f dx i r dx 2 K^B J . e (x 1 )A e+m (x 2 ) 

' ’ J ^,.1 CO J —CO 

77i:— 1 

A^(2xi + x 2 - 2x), (2.49) 

where 

K m _ e -4A(d+C 2 K?^3 5 (2.50) 

and (1.6.20) is used. With (2.7), (1.7.10) and (1.7.13) it is easy to show that 

?{ u $£,*>}( i ' k = °) = o- ( 2 - 51 ) 

The jump by the com P onen t) obtained using (2.5), (2.17), (2.45), (2.46), (1.7.31) 

and (1.7.41), can be written as 

f {u^ol^mr,} (*» = 0) = 4tt y dxi J dx 2 K^ c) A^t{x{)B m {x 2 )A* m _ i (x a + 2x 2 - 2x) 

+Z# c) +Z$?, (2.52) 

where 

A'{ 4c) = e - x [<Ci(Ci+&)*+2<|/ 3 ]^ 1 ^ 1 + £ a )f t} (2.53) 

77//^ = | 7r ^s / dxi f dx 2 f dx3C _ ^( 12 ^ +3 ^ + ^ +2l * 3 ^ 9 CiC3 1 + C3(2 £j - Ci)!,” 2 

J — OO J — OO J — DO 

B m {x 2 )A e . m (x 3 )A}_ :) (2x 2 + x 3 - 2x), (2.54) 
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tij$ = 8t rh, f dx, P dx 2 P dx 3 e~H^+^t+ 2 ^)/ 3 i k (;, [l - (3Ci + 2C a )€*er 2 ' 

J — OO J—OO J — oo *■ 

A J _^(x 2 ) J B m (a:3)A^_ / (x2 + 2x 3 - 2x), (2.55) 


& = 3C2 + ( 3 , = ii + 2Cl J = £7 — C2, 6 = 2 Cl + C21 (2.56) 


along with C 3 = x 2 — x 3 • Using (2.26) and integrating by parts, (2.54) and (2.55) can be 


rewritten as 


tip = 4vh s [ dx ,f dx 2 A^j C) B m (x 1 )A£_ TO (x 2 )i4^_ J (2x 1 + x 2 - 2x), (2.57) 

J —oo J — OO 

U { ni = 4jt h. f dx, P dx 2 \Kffl a B m (x,)A^(x 2 )A^_ e (2x, + x 2 - 2x) 

J—oo J — 00 L 

+ K n% A 3-i{ x \) B rn.{x2)A* m _ t {x 1 + 2x 2 - 2x)] , (2.58) 


where 


KP = ^(0, ( 2 ,( 3 ; 3, 1,1), Kpl = #(6,6,0; 1,1,1), Kpl = #(6,0,0; 1,|,2), 

(2.59) 

#(C., 6 , Cc; a, 6 , c) = ( 2 /d)e-UV 3 )[«?+{ 2 -(- 2 ) 4 }C?K fc 3 ] (2Ca + Cz) [( Cl /d)(3Ci + 26 C 2 ) 

+ jT Cl 3 exp {±(2aA/3) [(Cs =F 3Ci)Cl + 2 C 2 CI T (3/o)ftC*a] } {& - (2/a) ( 2 ^ + C 2 ) 

[l + A| a (3Ci + 2Cfe — C 3 + 4Cc)(2Ci — 2Ca + C& + C 3 )] }] , (2.60) 


and ia and |/ are defined in (2.32) and (2.56). 


3. Amplitude equations 


The non-equilibrium amplitude equations for the frequency-detuned resonant-triads 


are derived by substituting the velocity jumps across the critical layer obtained in §2 into 
the jump equations (1.7.55) and (1.7.56). 



By substituting (2.16), (2.18) - (2.20), (2.22), (2.36) and (2.40) along with (2.2) into 

(1.7.55) , we obtain the frequency-detuned amplitude equation for the oblique modes (see 
also (4.2)), for — J ^ j ^ J, 

( COS0 + ^) + {i “ Vo) - Kob] A)] = J_^dx l K x {x,x- i )B 1 { x x)A]{2x 1 - x) 

i J f s f Xl 

-—(sec 2 0) 2J / dxi dx 2 [K 2 (x,xi,x 2 )A. l (xi)Ae(x 2 )A{(xi + x 2 - x) 

4t _jJ— o° J — oo 

+ir 3 (*, *i, *2 )^(®iMj(®2MK x i + ^2 - *)] > (3.i) 

where x and X are defined by (1.12). The frequency-detuned plane-wave amplitude equation 
is similarly obtained by substituting (2.48), (2.49), (2.51) and (2.52) along with (2.2) into 

(1.7.56) , that is (see also (4.3)), for — J ^ j ^ J, 

4i v—' f s f Xi 

Bjx + {i Ox - Vo) ~ «2<i} Bj = 22 / dx ' dx2[K 4 (x,Xi,X2)Bj(xi)Ai(x2) 

T £_ j J — QO J-OO 

A}{ 2xi + x 2 - 2x) + Ks(x,Xi,X 2 )At(xi)B 3 (x 2 )A* t (xi + 2x 2 - 2x )] 

+(back — reaction — term). (3.2) 

The oblique-mode amplitude A 3 and the plane- wave amplitude B 3 are normalized by (1.13). 
At the leading order, r (= c/(r w Y c ) as defined by (1.7.18)) becomes unity and 9 given by 
(1.4) becomes tt/ 3. The transverse-coordinate origin shift r] 0 can be chosen to be zero as in 
(1.7.58). The analysis allows the value of J in (3.1) and (3.2) to be very large (as long as the 
magnitude of J\ in (1.6) is 0(1), see also §4). The upper and lower limits of the summations 
could also be arbitrary (i-e- £/=*)• The amplitude equations (3.1) and (3.2) (also (4.2) and 
(4.3)) are valid for any value of the local-growth-rate exponent r in the range 0 < r ^ 3. 
The analysis can be extended to include higher-order effects in the amplitude equations 
(3.1) and (3.2) with minor modifications in their coefficients, which may be required for 
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good comparison with experiments (see Wundrow, Hultgren & Goldstein 1994 for excellent 
comparison of the linear growth rate). 

The kernel functions, K\ for the parametric-resonance term, K 2 and K$ for the self- 
interaction term and and K 5 for the mutual-interaction term, are given by 


Ki = e~ 2X ( £ - Xl ) 3 / 3 (x-xi) 2 , (3.3) 

K 2 = Ji {3c) - 2r(sin 2 8 ) [hf/J - Af e > - A< 3j f) + hf^] , (3.4) 

A 3 = Kf d) - 2r(sin 2 8 ) [aJ5? - A'f e) - A< 3/) - K (3h) + A^J } ] , (3.5) 

K 4 = 2A< 4i > + r(sin 2 8 ) [k\] c) + K$ a ] , (3.6) 

A s = a{ 4c) + r(sin 2 8 )k\% (3.7) 


where A is defined by (1.12). The 0(1)- viscosity kernel functions for the self-interaction and 
mutual-interaction terms, which are obtained by substituting (2.21), (2.23), (2.28), (2.37), 
(2.38), (2.41) and (2.44) into (3.4) and (3.5) and by substituting (2.50), (2.53) and (2.59) 
into (3.6) and (3.7), respectively, are given Appendix A. The sum of K 2 in (A 1) and A3 in 
(A 2) is the same as the self-interaction kernel function of the single-frequency (or single- 
resonant-triad) interaction given in Wu et al. (1993). The last terms which involve C 1 C 2 in 
(Al) and (A 2) appear with opposite signs in K 2 and A3. 

The oblique modes react back on the plane waves and the corresponding back-reaction 
term is quartic in the oblique-mode amplitudes (Goldstein & Lee 1992; Wu 1992, 1995). 
The kernel functions for the back-reaction term will be analyzed in future. 

In the inviscid limit (A = 0), the kernel functions become 

Ki = (x-xj) 2 , (3.8) 
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K 2 = (l — 2r sin 2 o'j (x — Xi) 3 , (3-9) 

K 3 = (l - 2r sin 2 Oj (x - Xi)(x - x 2 ) [x - Xi + 2r(sin 2 0)(xi - x 2 )j , (3.10) 

Ki = 2 (l — r sin 2 (x — xi) 3 , (3-11) 

K$ = ^1 — t sin 2 9^ (x — xi)(x — x 2 )(2x — xi — x 2 ). (3.12) 

The linear growth rates of the oblique and plane waves, « 0 6 and n 2 d, are 


sec 6 


K 0 b — 


~a { m + £ (2a T KT w \) 


1/2 


, K2d = -affl + ! (° r *T w \) 1/2 , (3.13) 


r ■ 2 


4(1 + cos 2 6) [ r 

where the real constant k is a normalization parameter which can be chosen arbitrarily, 
r w is the total wall-shear stress (see (1.3.15) and (L3.16)), is defined by (1.12) and £ 
(= a/(r w c) as defined by (1.7.57)) is unity at the leading order. The last terms in (3.13) are 
due to the viscous Stokes-layer effect. Their magnitude is 0(<7 r / 2 ) when the critical layer is 
governed by the non-equilibrium dynamics, but it becomes 0(1) in the viscous limit where 
A = 0(c -r ). 

When the upstream flow is composed of a system of resonant-triads of linear instability 
waves (Goldstein & Lee 1992; Lee 1997a) the upstream boundary condition becomes, as in 
(1.7.59), 


A } -*• d^exp [(k 0 6 - ^./x)*] i Bj b : exp [(k 2cJ - ijx)^] as x -*• (3-14) 

where the complex initial amplitudes a 2 and b 3 are given by (1.7.60). In the later down- 
stream stage of the critical-layer evolution, the upstream boundary condition for the non- 
equilibrium amplitude equations in that local streamwise region is obtained by matching 
with the solutions of the preceding critical layer stage (Wundrow et al. 1994; Goldstein 
1994; Wu, Leib & Goldstein 1997; Lee 1997a). The amplitude equations (3.1) and (3.2) are 
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valid for any value of the local-growth-rate exponent r (it was assumed that 0 < r ^ 3 in 
Part 1), 

The stream wise evolution of the instability waves is determined by solving the frequency- 
detuned amplitude equations (3.1) and (3.2) along with the upstream conditions (3.14). The 
numerical solutions of these amplitude equations will be given in §9. The frequency-detuned 
self-interaction term is composed of two summations in this multi-resonant-triad-interaction 
case. The wavenumber differences between the instability waves of the Oth and jth resonant- 
triads (that are given by (1.5.21)) produce the terms that involve jx in (3.1) and (3.2). By 
putting J = 0, we can recover the previously obtained amplitude equations (Goldstein & 
Lee 1992; Wu 1992, 1995; Wu et al. 1993; Lee 1997a). The nonlinear part of the amplitude 
equations (3.1) and (3.2) includes only one mean-flow-dependent parameter A. 

In the later downstream stage of the non-equilibrium critical-layer evolution, the local- 
growth-rate exponent r becomes smaller than in the upstream non-equilibrium stage. The 
mean-boundary-layer flow develops on the long viscous length scale (x v given by (1.2. 1)). 
The boundary-layer thickness and the mean pressure gradient are nearly unchanged (within 
the order of approximation) over the region where the nonlinear interaction, or a series of 
nonlinear interactions, takes place. If we let the local-growth-rate exponent be r 0 in the 
(first) upstream non-equilibrium critical-layer stage, the Reynolds number, mean pressure 
gradient and frequency detuning are scaled by, from (1.8), (1.3) and (1.7), 

_L _ <7 3’-o+4 / \ 05 J! _ a T °~ l fi 0 , X = <r T °Xo, (3.15) 

where A e , £l 0 and \ 0 are 0(1) (of course they could be smaller). In the downstream region 
where the local value of r is smaller than r 0 (we have only used the first equation in (1.3) 
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for the simplicity), A, fi and x become 

A = <r 3 ( r °- r >A 0 , ji = a T °~ T ji 0 , X = ^°~ T Xc- (3.16) 

Therefore, the viscosity, mean pressure gradient and frequency detuning become less im- 
portant in the downstream non-equilibrium region where r < r Q (Wundrow et ai 1994). 

From (1-12), (3.13) and (3.16), we can show that the linear growth rates K 0 b and K 2 d 
become negligibly small when r is smaller than r 0 (where r Q ^ 3). Therefore, in this 
later non-equilibrium critical-layer stage, the instability wave amplitudes of the frequency- 
detuned resonant-triads are determined by the integro-differential equations (3.1) and (3.2) 
with the linear terms omitted and the kernel functions given by their inviscid limits. 


4. Large frequency-detuning limit 


If we put 


Aj = A,exp(!;x*)> B } = B : ex^(\j X x), 


(4.1) 


the non-equilibrium amplitude equations (3.1) and (3.2) can be rewritten as 


cos 


9 + { A J* ~ (i 7 ? 0 + K ° b ) A j} = ~ J dx 1 K 1 B J (x 1 )A*{2x 1 - x) 


1 o f ^ 

(sec 2 0) / dx i / dx 

4r ’ J - oo J-co 


(K 2 + K 3 )A } (xi)A :i (x2)A* j (xt 1 + X2-x)+ Ms 




,(4-2) 


Bfl ~ (i Vo + «2 d)Bj = / dx i / dx 2 K a B j {x x )A ] {x 2 )A :i {2xi + x 2 - 2x ) 

T J — OO J — OO L 

+ 1 ^1 — t sin 2 (sec 2 
f dx x f dx 2 f dx z \k 6 A :! (xi)A :i (x 2 )A :) (x 3 )A*(xi + x 2 + x 3 - 2x) + ^ Ad*] , (4.3) 

J — oo J — OO J — OO *■ J 


+K 5 A J (x x )B ] (x 2 )A*(x x + 2x 2 -2x) + ^ A4 m 
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where K\ - K$ are given in the previous section and Kq is given as, in the inviscid limit 
(see Wu 1995 for finite viscosity one), 

K$ = 2 r ^sin 2 9j £ 2 [^1 + 2r sin 2 o'j (tj + £ 2 ) 2 + (Cl + C 2 + C 3 ) 2 ] + ^1 — 2r sin 2 £1 
[(2Ci + C 2) 2 (3Ci + 2 C 2 + C 3 ) _ 2r (sin 2 {(1 + (Ci + C 2 ) (Ci + C 2 + C 3 ) (Ci ~ ( 3 )}] -(4-4) 

We have put £3 = x 2 — X 3 and £1 and {2 are defined in (2.3). 

The multi-mode-coupling terms M s and M m (-Mb is not derived in this paper), be- 
tween the yth and tth resonant-triads, are given by 

M s = e^ (} -^ g - Xl) K 2 A J (x 1 )A e (x 2 )At(x 1 + x 2 - x) 

+e$*b- e ><*-»)K s A t (x 1 )A J {x 2 )%(x 1 + x 2 - x), (4.5) 

M m = K 4 B ] (x 1 )A e (x 2 )Af(2xi + x 2 - 2x) 

+j*b-W*-**)KiA e {x 1 )Bj(x 2 )Afa 1 + 2x 2 - 2x), (4.6) 

When |x(j — t) | is large, we can show that, by integrating by parts in the inviscid case, 

f dx 1 r ixiM, - 96(1 ~ r &,|A,(*,)| 2 + 0(U(;-«r 5 ), (4.7) 

J-OO J-oo [xO “ VJ J-CO 

r dxi r dX 2 M m = r dx.iMx^+omj-^r 5 )- w 

J - 00 J-oc [x(j — t)j J-oo 

The effects of the multi-mode-coupling terms, M s and M m in (4.2) and (4.3), become neg- 
ligibly small when the frequency difference becomes very large (we assume similar behavior 
for the Mb term). 

In the real flow, there are infinite instabilities whose frequencies and wavenumbers 
are continuously varying. The individual critical layers are merged from the beginning 
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(of nonlinear process). The truncated finite summation in (1.10), (3.1), (3.2), (4.2) and 
(4.3) is an approximation to the infinite series. The instability waves whose wavenumbers, 
frequencies and linear growth rates can be scaled as in §1 are only included in the summation. 
The amplitudes of the other linear modes (i.e. of very high/low frequencies) have been 
assumed to be too small to play any major role in the nonlinear process. However, it is 
reasonable to expect that the growth of these other modes will be eventually governed by 
the frequency-detuned interaction in the downstream region where their nonlinear growth 
rates become much larger than the linear growth rates (that are different from (3.13)). 
The amplitude equations (4.2) and (4.3) (also (3.1) and (3.2)), that are obtained when the 
magnitude of Jx is of 0{a T ) (see (1.6)), can still be used, along with (4.7) and (4.8), for 
larger value of (but Jx still smaller than 0(1)). 

5* The singular solution of the amplitude equation 

The solutions of the amplitude equations (3.1) and (3.2) always develop a singularity 
at a finite downstream position x s as will be shown in §9 and was shown in Goldstein & 
Lee (1992), Wu (1992) and Lee (1997a) for the single- resonant-triad interaction. Near the 
singular point where the exponent r is smaller than in the upstream region, the linear terms 
and the viscosity effect become of higher order as explained in the previous section. The 
amplitudes are then determined by (3.1) and (3.2) with the inviscid kernel functions (3.8) - 
(3.12) and with the linear terms neglected. The solutions to these equations can be written 
as (Goldstein & Choi 1989; Goldstein & Lee 1992; Wu 1992), when x — ► 

A, = aj{x s - x) 3+i ^, Bj = bj/(x s - x) 4+2i ^ for 3 = J, (5.1) 
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where 'ipj is a real constant and a 3 and b 3 are complex constants. 

By substituting (5.1) into the inviscid amplitude equations (with the linear terms omit- 
ted), we can show that 

J 


( 1 \ 3 + iVh i b, i(l — 2r sin 2 9) la^l 2 _ , , . . 

( COS0 + JS?) = r^ Dp ^ j) ~ 4r(cos 2 6) ( } 

(5.3) 


4 + 2r0j_ 4i /t ; _2 \ a *\ 


' [ T' j ' 2 ~~ = (1 - t sin 2 6 ) ^2 fa) + (back - reaction - term), 

where the integrals D p , and D m are given in Appendix B. The above equations can be 

solved numerically to determine \a 3 |, \bj |, fa and the argument of 6,/a 2 . As in the single- 

resonant-triad case (Goldstein & Lee 1992), (5.2) and (5.3) fix only the argument difference 

arg(6^/a 2 ) and leave the individual arguments arg(a <? ) and arg(6j) undetermined. 

From (1.9) and (1.10), we can show that the streamwise velocities of the plane and 

oblique modes, in the inviscid wall layer, become of the order of the base mean flow in the 

downstream region where — x is 0{<J r ). 

From (5.1) along with (1.2) and (1.12), we can show that, as x -+ x s , 

A JX r +1 kr w aAjx r+1 kr w a (3 + i^) B JX TJrl kr w aB 3 x r +i (4 + 2i^,) 

, 2(x s -x) ’ 

(5.4) 


2 A 


2(x s -x) 5 Bj ~ 2 Bj 

Near the singular point, x s — x becomes small so that we can put 

x s - x - a n ( x s - x ) , 


(5.5) 


where the magnitude of x s — x is of 0(1). If the local growth rates are expressed as 
A 3X /Aj = 0(<7 r+1 ) and B JX /B 3 = 0(a r+1 ), then it follows from (5.4) and (5.5) that the 
local value f is given by 

r = r — n, (5.6) 
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which shows that the local-growth-rate exponent f becomes smaller as the distance x s — x 
becomes shorter. 

Since the critical-layer thickness in the non-equilibrium analysis is of the same order as 
the growth rate, the critical layer becomes thicker as the singular point is approached. The 
thicker critical layer leads to a wider range of instability-wave frequencies over which the 
frequency-detuned nonlinear interaction can occur (see (1.6)). All instability waves whose 
scaled frequency differences are nearly 0(1) can nonlinearly interact right before the scalings 
of this analysis break down, which occurs at the streamwise position where the critical layer 
is as thick as the inviscid wall layer and the unsealed distance from the singular point is of 
the order of the wavelength scale (Goldstein & Lee 1992; Wu et al. 1997). 

6. Viscous-limit quasi-equilibrium amplitude equation 

The quasi-equilibrium amplitude equations for the multi-resonant-triad interaction will 
be obtained by taking the viscous limit (Wu et al. 1993) of the finite-viscosity amplitude 
equations (3.1) and (3.2). Mankbadi, Wu & Lee (1993) and Wu (1993) derived the quasi- 
equilibrium amplitude equations for the single-resonant-triad interaction by considering the 
Tollmien-Schlichting waves in the upper-branch-scaling regime in the Blasius and favorable- 
pressure-gradient boundary layers, respectively. 

In the viscous limit as A — ► oo, the kernel K\ given by (3.3) becomes highly concen- 
trated around x = X\ and the parametric-resonance term N p in the oblique-mode amplitude 
equation (3.1) becomes (Goldstein & Lee 1993) 

ff, = ~B,(,i)A;(x). (6.1) 
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The viscous limit of the frequency-detuned self-interaction terms shows that the leading- 
order terms are produced by the second double-integral term (whose kernel function is K 3 ) 
in ( 3 . 1 ). Following Wu et al (1993), we can show that the viscous limit of the integral 
which involves the second term 2 h s f £ 2 d(s (. . .) in (A 2) becomes 

i / 1 \ 1/3 J rx 

~ ^ ( 18 ) r ^) J2 Mx) J^ d x 1 A. 1 ( x 1 )A.* t (x 1 ), (6.2) 

and the viscous limit of the integral which involves the fourth term C 2 : C 31 0? C 2 ) in 

(A 2 ) becomes 

ir / 1 \ 1/3 J rx 

^(tan 2 0 sin 2 0 ) (— ) r(|) £ M s ) J^dxiA^x^A^), (6.3) 

where T(l/3) = 2.6789 is the Gamma function. The above two terms in K$ are originated 
from Kj 3d \ Kj 3d \ and given by (2,23), (2.28), (2.41) and (2.44). The nonlinear 
interactions of the ^ 02 ^ an d ^ 02 components, as given by (1.7.26) and (1.7.30), 

are responsible for the velocity jumps of the (3d) and (3 h) components in (2.22) and (2.40) 
(Wu et al 1993). This is consistent with the results of the quasi-equilibrium critical-layer 
analyses by Mankbadi et al (1993) and Wu (1993) who showed that the spanwise-periodic 
mean-flow distortion is responsible for the velocity jump across the diffusion layer. The 
magnitude of the viscous limits of the integrals which involve the other terms of K 3 in 
(A 2 ) and the kernel K 2 in (Al) are all of 0(A“ 5 / 3 ) or smaller (Wu et al 1993). The 
frequency-detuned self-interaction term N s becomes 21 s A 00 , from ( 6 . 2 ) and (6.3), 

Ns = ~2X V3 ( 2 _2rsin20 ) ( tan2 0) (jg) r (|) M*) J_^dx 1 A ] {x l )A}{x l ). 

(6.4) 

Similarly, we can show that the mutual-interaction term N m in the plane- wave ampli- 
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tude equation (3.2) becomes in the viscous limit, 

N m = -^C m ^B } {x)\A e {x)\\ (6.5) 

where C m is a real constant that can be obtained by integrating the kernel functions K A 
and K 5 given by (A3) and (A 4) with A = 1 for the exponential terms and with that is 
redefined as not to include the last term A((i — ($)(. . .)(. . .) in (A 7), 

Cm=[ dx\ [ dx 2 fA r 4 (x,ii,x 2 |A = l;^*) + A' 5 (x,ii,i 2 |A = 1; J r 6 + )| . (6.6) 

J — oo J—c o J 

We can assume, from the back-reaction term in the single-resonant-triad case given by 
Wu (1995), that the viscous limit of the back-reaction term is 

N b = 0 (l/A N ) , (6.7) 

with N larger than 5/3, which has been indirectly proved by the quasi-equilibrium analyses 
of Mankbadi et al. (1993) and Wu (1993). 

With (6.1), (6.4), (6.5) and (6.7), the amplitude equations (3.1) and (3.2) become in 
the limit as A — ► oo, 

( COS * + cos 0/ Aj£ + ~ K ° b ) A l] = N P + ( 6 - 8 ) 

B : s + (i;x - «2 d) B } = 0, (6.9) 

where k ^ and K 2 d are given by (3.13) and we have put rj 0 = 0 (see (1.7.58)). 

If we rescale the viscosity parameter and the amplitudes as in Lee (1997a), 

A = a T A, Aj = a 2r/3 Aj, B 3 = <r T A } , (6.10) 
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the above equations (6.8) and (6.9) can be rewritten as (see also (6.17) and (6.18)), foT 


j = 


(cos 9 + [4,* + (i?X - Kob) ij = -^B 3 {x)A](x) 


; / N / i \ l/o w— , /*ar 

( ! " 2rsin20 ) ( tan 2 0 ) (jg) r(|) 53 ^(*) Ljx^xMX^ 

Bjx + (i;x - «2d) = o, 


( 6 . 11 ) 

( 6 . 12 ) 


with 

, k 2d = + | (kT w Xy /2 , (6.13) 

where the normalization parameter k can be chosen arbitrarily, r and f (as defined by 
(1.7,18) and (1.7.57)) are equal to one, \ an ^ m are & yen by (1*12) and r w is the total 
wall-shear stress (see (1.3.15) and (1.3.16)). The linear growth rates k 0 b and hid were 
obtained by substituting (6.10) into (3.13). As we expected, the viscous Stokes-layer effect 
becomes 0(1) in the viscous-limit linear growth rates. 

The streamwise evolution of a system of frequency-detuned resonant- triads that are 
composed of the Tollmien-Schlichting waves is determined by the quasi-equilibrium ampli- 
tude equations (6.11) and (6.12). The value of the exponent r is 3 for the Blasius boundary 
layer and 1 for the favorable- pressure-gradient boundary layer. These amplitude- equations 
(6.11) and (6.12) can also be directly derived from the frequency-detuned quasi-equilibrium 
critical-layer analysis (Mankbadi et al. 1993; Wu 1993). The quasi-equilibrium critical-layer 
scaling can be obtained by rescaling the non-equilibrium scaling given in §1 using (6.10) 
(see table 1 of Lee 1997a). The equations (6.11) and (6.12) (also (6.17) and (6.18)) for the 
frequency-detuned multi- resonant-triads reduce to the amplitude equations of Mankbadi et 
al (1993) and Wu (1993) if we put J — 0. 
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The plane waves grow linearly, i.e. 


B-, = b : exp [(k 2 d ~ ijx) *] , (6-14) 

in the entire quasi-equilibrium region. The streamwise evolution of A 3 can be obtained by 
solving (6.11) with the linear upstream condition 

Aj -*■ dj exp [(k o6 - %3X)x] • (6.15) 

The complex constants a 3 and b 3 are given by the equations similar to (L7.60). The results 
of the numerical computations are presented in §8. 

The growth of the frequency-detuned ToUmien-Schlichting waves in the upper-branch- 
scaling regime is governed by the quasi-equilibrium amplitude equations (6.11) and (6.12). 
However, this initial quasi-equilibrium stage will be eventually followed by the non-equilibrium 
critical-layer stage where the amplitudes are determined by the integro-differential equations 
(3.1) and (3.2) (Goldstein 1994; Wu et al. 1997; Lee 1997a). 

If we put 

Aj = ijexp(^xz), Bj = _Bjexp(yxx), (6.16) 

the equations (6.11) and (6.12) can be rewritten as 

(cos 6 + — ^) (i j5 - Kofclj) = -Xrb 3 {x)A 3 {x) 
-_ 7i (l-2rsin^)(tan 2 0) (-J T(§) A 3 {x) <fe 1 |A J (a: 1 )| 2 + Y M q ,(6.17) 

Bjx - = 0, (6.18) 

where 

M q = A e {x) r (6.19) 
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As in §4, we can show that 


Mq = ^—^ Aj(x)\A e (x)\ 2 + O([x0 - £)] 2 ), (6.20) 

in the large frequency-detuning limit (when \x(j — £)\ is large). 

7. Nonlinearly-generated low-frequency modes 

Goldstein & Choi (1989), Goldstein & Lee (1992) and Wu (1992) showed that the 
nonlinear interaction between oblique modes of the same frequency produces a spanwise- 
periodic mean-flow distortion. However, the nonlinear interaction between the frequency- 
detuned oblique modes generates low-frequency modes in addition to the spanwise-periodic 
mean-flow distortion. The frequencies of these nonlinearly-generated modes are equal to 
the frequency differences between the primary oblique modes, therefore, they are very low 
(of 0(cr r+2 )). From (2.8) and (L7.ll) along with (2.2), we can show that 

= [ U o%,e v dT ) = - 47rr (sin 2 6) [ dx x H x (x - x^A^x-^A^Xx), (7.1) 

AV o,2y,f(z) = / V Q%,t n d V = 4i7T r (sm 2 6) f dx x H 0 (x - x x )A J+ t(x x )A}(x x ), (7.2) 

J —CO 1 ’ ’ J — OO 

a <4, = n = o. (7.3) 

J — OO 

where 

H n {x)= ( X dxe~ 2l£3 ^{x-x) n , (7.4) 

Jo 

0, x and A are defined in (1.4) and (1.12), 77 is defined by (1.7.1) and r (given in (1.7.18)) is 
equal to one. We can also show from (7.1) and (7.2) that 

A Vjg/i) - -ijj [A£fg a , ( (i)] . ' (7.5) 


! i I 
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The streamwise evolution of A } is determined by solving the frequency-detuned ampli- 
tude equations (3.1) and (3.2), therefore, the velocity jumps (7.1) and (7.2) are determined 
from the nonlinear critical-layer interaction. In order to match with these nonlinearly gen- 
erated jumps, the low-frequency modes in the inviscid wall layer (where y = crY ) must be 
written as, for Y < Y c , 

U0 , 2 = a 3r+1 £ £ Re { ^ 2 ^(y, cos [2Z + (j + 2*)*i] , (7.6) 

2——2J i: — \ ' J 

2 J r _ (2)± _ 'l 

i>o ,2 = ^ r+3 £ £Re cos [2Z + (j + 2^i] , (7.7) 

J——2J t\—l ^ J 

along with similar forms for w 0 ,2 and p 0 , 2 , where ti(= a r+2 xst ) and Zi(= a T+1 \/3xaz/4) are 

defined by (1.5.24), x is defined by (1.12) and denotes 35 was S^ ven ™ 

(1.6.20). The spanwise velocity wo ,2 and the pressure po ,2 are continuous across the critical 

layer (see (7.3) and (1.7.16)). The magnitude of the streamwise component u 0 ,2 given by 

= - (2)± . . 

(7.6) is as big as that of the primary oblique mode (see (1.5.30) where Vq^-, 3,1 = U o, 2 ;.j,f eIJxr )- 
The nonlinearly-generated low-frequency modes are functions of the slow time t\ and the 
magnitude of their frequencies is of 0(a T+2 ). 

z ( 2 )± z ( 2 )± 

The shape functions U 0 2 ^ and V 0 2 . h( are discontinuous at the critical level Y = Y c . 
We can show from (1.6.5) - (1.6.7), (1.6.14) and (I.B2) that they have to satisfy the following 
equations 


= ( 2 )+ 


= ( 2 )- 


2k 3 . 


V'oY/Y - n + ) - U'oYAY = Y-) = U™,.,, 


M 


= ( 2 )+ 


( 2 )- 


k 4 ac . 


V'vY/Y = V« + ) - V'oYAY = Y~) = 


M 


(7.8) 

(7.9) 


where t and are given by (7.1) and (7.2), M is given by the phrase below 

(1.13) (also in (1.7.5)) and k is the normalization parameter introduced in (1.12). 
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In the upstream linear region where A } is given by (3.14), the critical-layer jump 


A UQ 2 - 3t e (7-1) becomes, as x — *■ — oo, 


2 &ob 2 JX 


(2A)V3 

( 2 Kob - Jjx) 


exp [(2 k 0 (, — %jx) x , (7.10) 




where Hi is the Airy function (Abramowitz & Stegun 1965, p.448). In the inviscid limit 
where A = 0, (7.10) can be written as 

= -47rr(sin 2 9) (2n ob - %]x) 3 exp [(2 k o6 - £jx)x] . (7.11) 

The transverse velocity jump AVqj-jj can be obtained by differentiating (7.10) and (7.11) 
with respect to x as shown by (7.5). The streamwise growth rates of A and A V$ e 
in the linear region are twice the linear growth rate of the oblique mode. 

Near the singular point the viscous effect becomes small and the amplitude is given by 
(5.1). We can show that the jump (7.1) becomes, when x — ► 

au (2) 47rr(sin 2 fl)a J - K a; 1 19 > 

°’ 2y ’ < (3 + iA^)(4 + iA^)(5 + iA^) (x s - x) 3+iA ^' ’ ( J 

where we have put A = ^j+t — 

The numerical evaluation of the velocity jumps A U^jt anc ^ AV ^ ^ will be presented 
in §9 along with the numerical solutions of the amplitude equations. 

The analysis in the critical layer only determines the critical-layer jumps of the stream- 
wise and transverse velocities of the low-frequency (spanwise-periodic) modes. The com- 
plete solutions of the nonlinearly-generated low-frequency modes can be obtained from the 
multi-layer analysis (Wu 1993). When r = 3 and A defined by (1.8) is 0(1) in the single- 
resonant-triad-interaction case (Goldstein & Lee 1992), we need to consider the steady 
spanwise-periodic (mean-flow-distortion) mode in the potential region where y = y/cr, the 


32 


main boundary layer where y = 0(1) and the viscous wall layer where y = a s Y in addi- 
tion to the inviscid wall layer of 0(a) and the critical layer of 0(cr 4 ). The viscous wall 
layer, which is thicker than the viscous Stokes layer for the oblique and plane waves of the 
resonant- triad, is required in order to satisfy the wall boundary conditions. (The inviscid- 
wall-layer solutions of the streamwise and spanwise components become singular at the 
wall.) The three-dimensional boundary-layer equations in the viscous wall layer must be 
solved with the no-slip boundary conditions at the wall and appropriate boundary condi- 
tions on the upper edge of the layer in order to match with the solutions in the inviscid 
wall layer. The transverse velocity on the upper edge of the viscous wall layer is determined 
internally by the boundary-layer equations themselves and so its value can not be given as 
a boundary condition as was done in Wu (1993). It is interesting to note that the stream- 
wise and transverse velocity components in the inviscid wall layer, at leading order, become 
equal to zero where Y > Y c but non-zero where Y < Y c (Lee 1997c). Thus, the nonlinear 
interaction between oblique modes of the same frequency produces, at the leading order, the 
spanwise-periodic mean-flow distortion only below the critical layer. The non-zero velocities 
at y = y c _ will be determined by matching with the critical-layer jumps. 

Since the critical-layer jump of the spanwise-periodic mean-flow distortion is deter- 
mined from the solutions of the amplitude equations, which become singular at a finite 
downstream position, the nonlinearly-generated spanwise-periodic mode becomes very large 
near the singular point. When the value of the local-growth-rate exponent r becomes 9/10, 
the streamwise velocity of the mean-flow distortion becomes as big as the base mean flow 
near the wall. Therefore, the nonlinearly-generated spanwise-periodic mode and the base 
mean flow start to interact nonlinearly in the nonlinear viscous wall layer whose thickness 
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is of 0(cr 37 / 10 ) (Lee 1997c). This nonlinear viscous wall layer will be separated into two 
layers (inviscid nonlinear layer and the viscous nonlinear layer) in the later downstream 
region. The details of the evolution of this nonlinearly-generated mode will be presented in 
a forthcoming paper. 

8. Numerical solutions of the quasi-equilibrium amplitude equations 

The quasi-equilibrium amplitude equations derived in §6 are solved numerically. The 
frequency-detuned amplitude equations (6.11) and (6.12) determine the streamwise evolu- 
tion of a system of frequency-detuned resonant-triads of the Tollmien-Schlichting waves. 

A predictor-corrector method (see Gear 1971) is used to solve the oblique-mode ampli- 
tude equation (6.11) subject to the upstream condition (6.15), along with the linear plane 
wave given by (6.14). As in Lee (1997a), the Adams-Bashforth method was used for the 
predictor step and the Adams-Moulton method was used for the corrector step. The typical 
streamwise grid size is 0.002 and the streamwise integration starts about x = —20. 

The numerical results in this section are obtained when A = 1, k 0 b = 4/5 and k^d = 1. 
The exact linear growth rates can be obtained from (6.13) for the specific problem. We 
have also put Q — i r/3 and r = 1 (see (1-4) and (1.7.18)). The initial amplitudes of the 
plane and oblique modes are b : — 1 and a 3 = 0.01 exp(i7r/4) for all j. We have put the 
argument of the initial oblique-mode amplitude to be 7r/4 radian since it gives the most 
effective parametric-resonance growth of the oblique modes as shown in Lee (1997a). 

The general behavior of the solutions of the quasi-equilibrium amplitude equations in 
this multi-resonant-triad case is the same as in the single-resonant-triad case (Mankbadi 
et al. 1993; Wu 1993; Lee 1997a). The upstream linear growth of the oblique modes is 
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enhanced due to the parametric-resonance effect when the plane wave amplitudes become 
sufficiently large. The self-interaction between the frequency-detuned oblique modes causes 
the oblique-mode amplitudes to become saturated, which then grow linearly with oscilla- 
tions. Meanwhile, the plane wave continues to grow linearly in the entire region. 

Figure 1 shows the results of the nonlinear interaction between two frequency-detuned 
resonant-triads when j = 0 and 1. The solution of the single resonant-triad interaction is 
plotted as the dotted curve. The linear oblique-mode amplitude is plotted as the dot-dashed 
curve. Both plane-wave amplitudes Bo and B\ grow linearly as given by (6.14). In figure 
la, the result when the frequency-detuning factor x is equal to zero is plotted as the dashed 
curve (Ao and A\ are identical). The saturation amplitude is smaller than that of the single 
resonant-triad interaction because of the doubled self-interaction effect when x = 0* F° r 
the non-zero values of x, the saturation amplitude of Ao is larger than that of A\. The 
frequency of the 0th resonant-triad is lower than that of the 1st resonant-triad (see (1.6)). 
When the frequency detuning is very large the multi-mode-coupling term M q in (6.17) 
becomes very small. Thus, the resonant-triads grow almost independently of each other 
and the amplitudes Ao and A\ become close to that of the single resonant-triad interaction 
(within the range of the computation) as shown in figure Id. 

The numerical results of the three and five resonant-triad interactions are given in 
figures 2 and 3, respectively. The straight lines for the plane wave amplitudes are not 
plotted. The oblique-mode amplitude for the single resonant-triad interaction is plotted as 
the dotted curve. Both figures 2 and 3 (and figure 1) show that the saturation amplitude 
of the oblique mode of lower frequency (i.e. the one that has the lower value of j) is larger 
than that of higher frequency. 
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Since the effect of the spanwise-wavenumber (or propagation angle) detuning is not 
included in this analysis, the oblique-mode amplitudes |A,| are the same for all j in the 
parametric-resonance region. The spanwise-wavenumber detuning can be included by re- 
laxing the condition (1.5.16). The instability waves whose (scaled) spanwise-wavenumbers 
are detuned by 0(a T ) share the same critical layer and, thus, nonlinearly interact. 

9. Numerical solutions of the non-equilibrium oblique-mode amplitude equa- 
tions 


The results of the numerical computations of the frequency-detuned oblique-mode am- 
plitude equation (3.1) will be presented in this section. The critical-layer velocity jumps 
(7.1) and (7.2) of the nonlinearly-generated low-frequency (spanwise-periodic) modes will 
also be evaluated. 

We have observed in the single-resonant-triad-interaction analyses by Goldstein & Lee 
(1992) and Lee (1997a) that the self-interaction between oblique modes is mostly responsible 
for the explosive growth of the amplitudes near the singular point. Therefore, we will 
investigate the effect of the frequency- detuned self-interaction term first. There are, of 
course, many important shear flows (i.e. compressible free shear layer, supersonic boundary 
layer and others) where the the oblique mode is the most unstable wave and the self- 
interaction is the dominant nonlinear interaction. 

It is convenient to choose the normalization parameter k, which was introduced in 
(1.12) (or (1.7.1)), to be 


K = 


2 \tYcO_ 
T?,C 


z (2) 

1 1M> 


(9.1) 
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and we have shown in (1.4) and (1.7.18) that 


0 — 7r/3, r = 1. (9.2) 

If we put (1.7.58), (1.12), (9.1) and (9.2) along with B 3 = 0 into (3.1) and (3.13), the 
frequency-detuned non-equilibrium amplitude equation of the oblique modes can be rewrit- 
ten as, for — J ^ j ^ J, 

A 3 x = (§ - \]x) A 3 - |i f dx i/ dx 2 [K 2 {x,x- l ,X 2 )Aj{x 1 )Ae(x 2 )A* e (x 1 + x 2 - x) 

' 1-ZZ.—J J ~°° J —OO 

+K 3 (x, zuZ 2 )At(xi)A 3 (x 2 )Ai(xi + x 2 - x )] , (9.3) 

where the kernel functions iv 2 and #3 are given by (3.9) and (3.10) in the inviscid limit. 
The above equation (9.3) will be solved with the linear upstream condition (3.14), which 
can be rewritten as 

Aj — ► djexp [(| — as £ — » — 00 . (9.4) 

The effects of the magnitude of the initial amplitudes \a 3 \ and the frequency-detuning 
factor x on the streamwise evolution of the amplitudes A 0 will be shown in the following 
subsections when A = 0. 

As in before (Goldstein & Lee 1992; Lee 1997a), the Adams-Moulton method (see 
Gear 1971) is used to advance the solutions downstream from the prescribed upstream 
linear state. The double-integral term on the right-hand side of (9.3) is computed using the 
Newton and Cotes’ integration formula (see Kopal 1961). The numerical computation starts 
about x = -15 and the typical streamwise grid size is 0.01 (or 0.04 in the 51-pair-interaction 
case). 
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9.1. Inviscid two-, three- and five-pair interactions 


The numerical solutions of (9.3) with (9.4) in the inviscid case (A = 0) are presented when 
the instability waves are composed of the frequency-detuned two, three and five pairs of 
oblique modes. 

Figure 4 shows the results of nonlinear interaction between two pairs of oblique modes, 
i.e. j — 0 and 1, when the initial amplitudes are the same do = d\ = 1. The solution of the 
single-pair interaction is also plotted as the dotted curve. When the frequency-detuning 
factor x is equal to zero, the 0th and 1st amplitudes Ao and A\ become identical as shown 
in figure 4a. They become singular at earlier streamwise position than in the single-pair- 
interaction case because of the doubled self-interaction effect. This solution when x = 0 is 
replotted as the dashed curve in figures 4b, 4c and 4d. 

For the nonzero values of x, the initially same linear amplitudes Ao and A\ start to 
diverge when the effects of the frequency-detuned self-interaction become large as shown 
in figures 4b and 4c for x ~ 2 and 5, respectively. The multi-mode-coupling term M s in 
(4.2) becomes 0(1) and the frequency-detuned oblique pairs start to interact nonlinearly 
when the local-growth- rate parameter o T becomes as large as the scaled Strouhal number 
difference S\ — s (that is equal to (sj — s)jo 2 , see (1.5)). When x = 2 shown in figure 4b, 
the amplitudes Ao and A\ become nonlinearly interactive at earlier streamwise position (or 
at smaller growth rate) compared to the x = 5 case (figure 4c). In figure 4c, the growth 
rate of \Ao\ in the streamwise region where 1 < x < 3.3 is enhanced mainly due to the 
self-interaction between the oblique modes of the Oth pair. Similarly, the self-interaction 
between the oblique modes of the 1st pair is responsible for the enhanced growth of \A\\ 



there. Since the frequency detuning is relatively large (x — 5), the coupled interaction 
between two pairs occurs in the later downstream region, x > 3.3, where the growth rates 
become sufficiently large. Figures 4b and 4c show that \Aq\ is larger than \Ai\ (except 
1 < x < 3.3 in figure 4c). The frequency of the Oth pair is lower than that of the 1st pair 
as given by (1.6). 

When the frequency difference is very large, the coupled interaction between two oblique 
pairs does not occur until the growth rates also become very large. The results of x = 100, 
given in figure 4d, show that |Ao| and \Ai\ grow almost identically as in the single-pair- 
interaction case. The solid curve is not distinguishable from the dotted curve that is the 
solution of the single-pair interaction. 

It can be observed in figure 4 that both amplitudes become singular at the same 
streamwise position for all values of x. The frequency difference between two pairs delays 
the singularity. 

The streamwise evolution of the oblique-mode amplitudes are plotted in figure 5 for 
the three-pair interaction when j = 0, 1 and 2. The results are independent of the specific 
value of j of the center frequency, thus, j = — 1, 0 and 1 instead of the current choice 
will produce the same results as those in this figure. The initial amplitudes are the same, 
do = a\ = 0,2 = 1, as in the two-pair-interaction case given in figure 4. The solution of the 
single-pair interaction is plotted as the dotted curve and that of the three-pair interaction 
with x = 0 is plotted as the dashed curve. Similar to the results in figure 4, the coupled 
interaction between frequency-detuned pairs occurs at smaller amplitude when x = 2 (figure 
5a) than in the x = 5 case (figure 5b). When x = 5 the singularity occurs at later streamwise 
position compared to the case when x = 2. The magnitude of the Oth amplitude \Aq\ (of the 
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lowest frequency) is the largest and that of the 2nd amplitude \A 2 \ (of the highest frequency) 
is the smallest in the frequency-detuned nonlinear-interaction region. The amplitude of the 
lowest frequency, Aq, is likely to grow monotonically. The growth rate of \A 2 \ in figure 
5a becomes negative in a short period. It is interesting to observe in figure 5b that |Ao| is 
almost identical to the single-pair solution, but \A \ | and \ A 2 \ are smaller than that. In figure 
5b, the singularity occurs at nearly the same position as that of the single-pair solution. 

Figure 6 shows the results of the five-pair interaction when do = d\ — a 2 = = a 4 = 1 

with the results of the single-pair interaction (dotted curve) and the five-pair interaction 
when x = 0 (dashed curve). The amplitude of the Oth pair |j4o| is larger than the others, 
however, the highest-frequency amplitude \A 4 \ is not necessarily the smallest one (cf. figure 
5). The amplitude \A 4 \ still belongs to the group of smaller amplitude in the frequency- 
detuned-interaction region. Note that when \ = 4, shown in figure 6b, all amplitudes are 
smaller than the amplitude of the single-pair interaction and the singularity occurs even 
at later streamwise position than in the single-pair case. The singular point moves closer 
to that of the single-pair interaction for larger values of \ (the results are not shown). 
The solutions of ten-pair interaction show the similar behavior although the results are not 
presented. 

Figure 7 shows the solutions of the two-pair interaction when the initial amplitude of 
the lower frequency (fi 0 — 1) is larger than that of the higher frequency (a x = 0.8 and 
0.5 for figures 7a and 7b, respectively) and \ — 2. The single-pair solutions are plotted 
as the dotted curves for the respective initial amplitudes. (The dotted curve that has the 
same initial value as the curve (i) in figure 7a, for example, is the single-pair solution with 
do = 1 and the other dotted curve is the single-pair solution with the initial amplitude of 
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0.8.) The frequency-detuned interaction between two pairs causes the growth rate of the 
higher-frequency pair to become negative in a short period in both cases. The magnitude 
of the 0th amplitude \Aq\ is always larger than \A\\ 9 but both amplitudes become singular 
at the same streamwise position. 

The oblique-mode amplitudes are plotted in figure 8 when di(= 1) is larger than a 0 , 
which is equal to 0.8 or 0.5, and x = 2. Since the frequency-detuned self-interaction between 
two pairs initially enhances the growth of the lower-frequency amplitude Ao and reduces 
that of the higher- frequency amplitude A\ (as in figure 4), \Aq\ becomes larger than \A\\ in 
the later downstream region when do = 0.8 and d\ = 1 as shown in figure 8a. 

The amplitudes Ao, A\ and A2 of the three-pair interaction plotted in figure 9 show 
that |Ao| is always larger than | A2I although they have the same initial value (fio = 62 = 0.7 
and 0.5; d\ = 1). Figure 9a shows that the initially small 0th amplitude J Ao| becomes larger 
than the 1st amplitude \Ai\ in the downstream region similar to the two-pair-interaction 
solutions given in figure 8a. 

Figure 10 shows the effect of x on the streamwise evolution of the amplitudes in the 
three-pair-interaction case when ao = a2 = 0.7 and d\ = 1. The single-pair-interaction 
solutions with the respective initial amplitudes are plotted as the dotted curves. The onset 
of the singularity is delayed as x is increased from 0 to 3 as shown in figure 10a. When 
X = 5, the singularity occurs at earlier streamwise position compared to the case when 
X = 3. The amplitude ratios |Ao/Ai| and |A 2 /Ai| remained the same when x = 0. Figure 
lOe shows that the amplitude ratio of the higher-frequency mode IA2/A1I is smaller in the 
downstream region than the upstream value of 0.7 for all x considered. The ratio |Ao/Ai| 
for the lower-frequency mode becomes larger than the upstream value when x is equal 
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to 1, 2 or 3, but it becomes smaller than 0.7 when x = 5 as shown in figure 10c. The 
amplitude \A 2 \ of the higher-frequency pair exhibits more oscillatory behavior than \Aq\ of 
the lower-frequency pair as given in figures 10b and lOd. 

9.2. Inviscid 51-pair interaction 

The system of 51 frequency-detuned amplitude equations (9.3) where J = 25, i.e. j = 
—25, . . . , 0, . . . , 25, are numerically solved when A = 0. As we mentioned before the specific 
value of j of the center pair (0 in this case) does not affect the results. The initial amplitudes 
dj in (9.4) is given by the Gaussian function as 

a : = ^ exp - (; - n a ) 2 /50] + ^ for - 25 $ ; ^ 25, (9.5) 

in the single-peak cases (figures 11 and 12). The second Gaussian function 

exp [- (j - n 2 f /50] + (9.6) 

is superimposed on top of (9.5) for the initial amplitudes in the double-peak cases presented 
in figures 13 to 15. 

The Strouhal number difference between the jth and the Oth pairs of oblique modes is, 
from (1.6) and (1.12), 

Sj-s = a r ]XKT w s, (9.7) 

where k is given by (9.1). 

Figure 11 shows the numerical solutions when x is 0.2 and ni in (9.5) is zero. The mag- 
nitude of the amplitudes |Aj| versus j at different stream wise positions are plotted in figure 
11a. The streamwise evolution of the magnitude of the -20th, -10th, Oth, 10th and 20th 
amplitudes, |A(_ 2 o)|, |^4(-io)li |zlio| and |A 2 o|, respectively, is plotted in figure lib. 

! I!: 
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Figures 11a and lib show that the amplitudes of the instability waves, especially those of 
the lower- frequency modes (i.e. with smaller 7), are greatly enhanced due to the frequency- 
detuned self-interaction. The initially very small amplitudes of the -25th to -15th pairs 
become large as the instability waves evolve downstream. The ratio |A(_ 15 )/A(_ 25 )|, which 
is nearly equal to one in the upstream region, becomes large in the downstream region. The 
initial symmetric shape of the Gaussian function becomes asymmetric and the maximum 
amplitude occurs at lower frequency (at 7 = — 2 when x ^ 1) in the downstream region. 
In the upstream positions, up to about x = 0.5, the magnitude of the -20th amplitude is 
nearly equal to that of the 20th pair as shown in figure 11b. However, at later downstream 
positions (for example, at x = 2.52 and 2.96), |A(_ 2 o)l becomes much larger than |A 2 o| as 
well as | A(_ 10 )| becomes larger than j A10! as given in figure 11a. The amplitude | A20I still 
larger than the linear one at these two downstream positions (figure lib). All amplitudes 
become singular at the same finite downstream position. 

The effect of larger frequency-detuning factor is shown in figure 12 when x — 0.5 and 
7ii = 0. Because of the larger detuning, the frequency-detuned nonlinear interaction occurs 
at later downstream position in figure 12b compared to the previous x = 0.2 case given in 
figure lib. Figure 12a show r s that the frequency-detuned interactions enhance the growth 
of the lower-frequency pairs more than that of the higher-frequency pairs similar to the 
previous results. The frequency where the maximum amplitude occurs becomes lower as 
the oblique modes evolve downstream. The amplitudes become singular at later downstream 
position compared to the smaller frequency-detuning-factor case given in figure 11. 

Figure 13 shows the evolution of the double-peak initial amplitudes when x — 0.2, 
n\ = 5, n 2 = —5 and C 2 = 1. The amplitude |A(_ 5 )| is equal to | A5 1 in the upstream 
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region, but the former becomes larger than the latter in the downstream region. The 
amplitude of the -6th pair is the largest at x = 2.8. Figure 14 shows the results when the 
higher- frequency maximum of the initial amplitude, a$, is larger than the lower-frequency 
one, d 0 , and x = 0*2, n\ — 5, — 0 and C 2 = 0.8. Although the 0th amplitude is smaller 

than the 5th amplitude in the upstream region, |A 0 | eventually becomes larger than \A$\ 
in the later stage as was also observed in figure 8 in the two-pair interaction. Figure 15, 
which are the numerical results when x — 0.2, n\ = 0, = 5 and C 2 = 0.8, shows that 

the initially small higher-frequency peak (at j = 5) becomes less and less noticeable as 
the instability modes evolve downstream. The results in figures 13 to 15 also exhibit the 
enhanced growth of the lower-frequency pairs. 

The results of these frequency- detuned multi-pair interactions show that all instability 
modes become singular at the same finite downstream position. However, as we all know, 
the flow does not support extremely large instability waves. Near the singular point (where 
x s — x = 0(<j - 1 )), the amplitudes of the instability modes become as large as the base 
mean flow in the inviscid wall layer (see (1.5.30) and (5.1)). The next stage will be the 
triple-deck stage a s shown by Goldstein & Lee (1992) and Wu et al (1997). The flow in 
this stage is governed by the equations which are elliptic in the streamwise direction. The 
flow may enter the triple-deck stage before we can fully observe the explosive growth of the 
instability modes predicted by the frequency-detuned non-equilibrium amplitude equations 
of this analysis. The explosive growth occurs in a very short streamwise distance as the 


numerical results indicate. 



9.3. Critical-layer jumps of the nonlinearly-generated low-frequency modes 

When A = 0, the velocity jumps across the critical layer of the low-frequency modes, (7.1) 
and (7.2) can be rewritten as 

AU o%,t = “§* [ dxi(z - Xi) 2 A, + *(xi)A;(xi), (9.8) 

AV o%,t = 3i7r / dx x {x - Xi)A 3 +t{xi)A* t (x x ), (9.9) 

where we have used (9.2). The jumps and AVq^-j* are obtained by integrating 

(9.8) and (9.9) using the numerically calculated oblique-mode amplitudes A r 

Figure 16 shows the numerical results of the single-pair interaction when do = 1* The 
nonlinear interaction between a pair of oblique modes of the same frequency produces the 
steady spanwise-periodic mode. The oblique mode amplitude Aq that was plotted as the 
dotted curve in figure 4 is replotted as the same dotted curve. The streamwise growth rates 
of AU o%o,o and AV™ o are twice that of the oblique mode in the linear upstream region as 
given by (7.11). Both streamwise and transverse jumps of the spanwise-periodic mean-flow 
distortion become singular at the same singular point of the oblique-mode amplitude. 

The nonlinear interaction in the critical layer between three frequency- detuned pairs 
of oblique modes produces nine low-frequency (spanwise-periodic) modes. The numerical 
results of the three-pair interaction when do = di = d 2 = 1 are given in figure 17 for x — 2 
and in figure 18 for x = 5. The streamwise evolution of the amplitudes was plotted in figure 
5. If we put the subscripts j of A 3 to be 0, 1 and 2, we can show from (7.6) and (7.7) that 
the subscripts (j,£) of the non-zero AU^yt and AV$ t are (-2,2), (-1,1), (-1,2), (0,0), 
(0,1), (0,2), (1,0), (1,1) and (2,0). Figures 17 and 18 show that |A{7 q 2 2 ^| and [AT^-j*! 
with larger value of |j| are smaller than those with smaller |j| in the linear upstream region 
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as was predicted by (7.11). The mode with j = i — 0 is steady and its critical-layer jump is 
larger than the other components in most of the region. As in figure 5 the results in these 
figures are independent of the specific value of j of the center frequency. 

10. Concluding remarks 

A system of resonant-triads can interact nonlinearly between themselves in the common 
critical layer if their frequencies (of the fundamental plane waves) are different by a factor 
whose magnitude is of the order of the growth rate multiplied by the wavenumber of the 
instability waves. The long- wavelength small-growth-rate instability modes in boundary 
layers with and without mean pressure gradient are analyzed using the generalized scaling 
of Lee (1997a). 

In this part of the study, the system of partial differential critical-layer equations along 
with the jump equations given in §7 of Part 1 is solved analytically to obtain the frequency- 
detuned amplitude equations. The amplitude equations are similar to those obtained by 
Goldstein & Lee (1992) and Wu (1992) for the single resonant-triad interaction. However, 
in the multi-resonant-triad-interaction case, the frequency- detuned self-interaction term in 
the oblique-mode amplitude equation and the frequency- detuned mutual-interaction term 
in the plane- wave amplitude equation are composed of 2(2J + 1) terms (see (3.1) and (3.2)). 
The kernel functions for the self-interaction terms are divided into two parts. 

When the scaled- Strouhal- number difference between resonant-triads is larger than the 
local-growth-rate parameter <7 r , the multi-mode-coupling terms M s , M m and Mb in (4.2) 
and (4.3) become negligibly small (see (4.7) and (4.8)) and the instability waves of an 
individual resonant-triad grow independently of the other resonant- triads. 


The effect of the self-interaction between frequency-detuned pairs of oblique modes 
was investigated by solving the non-equilibrium oblique-mode amplitude equations in the 
inviscid limit. As given in the previous section, the growth of the lower-frequency oblique 
mode is more enhanced than that of the higher-frequency in most cases. It is shown that 
all amplitudes become singular at the same finite downstream position regardless of their 
initial values. The frequency-detuning effect delays the occurrence of the singularity. The 
singularity in an appropriately frequency-detuned multi-pair interaction occurs even at later 
downstream position than in the single-pair interaction although the sum of the initial 
amplitudes in the former case is much larger than in the latter case. 

Though the analysis was carried out with the non-equilibrium critical-layer scalings, 
the frequency-detuned quasi-equilibrium amplitude equations can be obtained by taking 
the viscous limit of the 0(l)-viscosity non-equilibrium amplitude equations. The numerical 
computation of the quasi-equilibrium equations shows that the saturation amplitude of the 
oblique mode of lower frequency is larger than that of higher frequency. The frequency- 
detuned resonant-triads of the Tollmien-Schlichting waves in the initial critical-layer stage 
is governed by the quasi-equilibrium amplitude equations, however, the later downstream 
stage will be eventually governed by the non-equilibrium dynamics as shown by Goldstein 
(1994), Wu et al. (1997) and Lee (1997a). 

Corke & Gruber (1996) and Liu & Maslowe (1998) show that the resonant-triad theory 
of Goldstein & Lee (1992) is in good agreement with their experimental and numerical 
results in adverse-pressure- gradient boundary layers. Corke & Gruber (1996) show that 
faster linear growth of the plane wave accelerates the parametric- resonance growth of the 
subharmonic oblique modes and induces larger saturation amplitudes of the oblique modes 


47 



compared to the Blasius case given in Corke k Mangano (1989). Using the parameters 
in Corke k Gruber (1996), we can show that the parametric-resonance effect is already 
much larger than the linear-growth effect (in the oblique-mode amplitude equation) at their 
first measured streamwise position. Faster (than linear) growth of the oblique modes at the 
beginning was observed in the experiments and also confirmed by the numerical simulations 
of Liu k Maslowe (1998). 

The numerical solutions of the non-equilibrium amplitude equations show that the ex- 
plosive growth of the instabilities occurs in a very short streamwise period near the singular 
point. Therefore, in most flow conditions for the resonant-triad interaction, the explosive 
growth may not be clearly distinguishable before the scalings of this analysis break down 
as the experimental (Corke k Gruber 1996) and numerical (Liu k Maslowe 1998) results 
indicate. However, the enhanced growth of the plane wave, especially in Figure 10b of Corke 
k Gruber (1996), proves that the oblique modes are large enough to produce the mutual- 
interaction and back-reaction effects on the plane wave, thus, there exists self-interaction 
effect. The explosive growth is the downstream asymptotic behavior of the self-interaction 
(or fully-coupled interaction) effect. In fact, as the finite- viscosity numerical solutions of Lee 
(1997a) (for example, his Figures 3 and 4) show, the fully-coupled interaction initially, in a 
very short period, reduces the growth rates of the oblique modes, which is more apparent 
when the viscous effect is large. The self-interaction effect can also be identified by the 
spanwise-periodic wall shear stress and the generation of the spanwise-periodic mean-flow 
distortion and nonlinear (very) low-frequency modes. The saturation of amplitudes, shown 
both experimentally (Corke k Gruber 1996) and numerically (Liu k Maslowe 1998), will 
be governed by subsequent nonlinear dynamics but the frequency-detuned self-interaction 
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between saturating/saturated oblique modes and neighboring frequencies (of small ampli- 
tudes) may still be active there to accelerate the growth of the latter. 

The present analysis is more focused on the frequency-detuned self-interaction and the 
summation does not appear in the parametric-resonance term in the oblique-mode ampli- 
tude equation (3.1). We can extend the analysis, by relaxing the condition (1.5.16), to 
include the spanwise- wavenumber detuning. It can be shown that the parametric-resonance 
interaction between a band of fundamental plane waves (w ± Aw/) and a band of subhar- 
monic oblique modes (u;/2 ± Ao; s ) can accelerate the growth of the oblique modes whose 
frequencies are within u>/2± (A w/ + Au? s ). The growth of a wide frequency range of oblique 
modes will be enhanced by subsequent interactions, for example, between the plane waves 
and the oblique modes of uj 2 ± (Au;/ + Au s ). Kachanov & Levchenko (1984) (also in the 
review by Kachanov 1994) showed that the (subharmonic) resonance is very wide in the 
frequency spectrum and it can amplify even (quasi-) subharmonic modes whose frequency 
detuning is close to one half the subharmonic frequency. 

In an incompressible boundary layer, the plane wave is the most unstable mode. How- 
ever, three-dimensional modes of small propagation angles are more unstable than the plane 
wave off the most unstable frequency (see (3.13)). Therefore, in a natural flow, the linear 
upstream flow may be dominated by an oblique mode whose propagation angle is relatively 
small. The resonant-triad analysis of this paper (also of Goldstein h Lee 1992, Wu 1992 
and Lee 1997a) can be generalized for three (pairs of) oblique modes instead of the usual 
fundamental plane wave and a pair of subharmonic oblique modes. The three (pairs of) 
oblique modes, whose streamwise and spanwise wavenumbers are a t j 2 and /?,, respectively, 
share the same critical layer, provided (aj 2) 2 + /? 2 are the same for all i — 1, 2 and 
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3. If we put the propagation angle of the first mode to be 9\ (= tan” 1 (2/3 1 /a 1 )), where 
0 $ 0i < 7r/6, the second and third ones must be 0 2 = 7r/3 — 9\ and 0 3 = 7r/3 + 9\ in order 
for the resonance-interaction to occur. When Q\ is equal to zero, 62 and 63 become 7r/3 and 
the three-oblique- pair resonant-triad reduces to the traditional resonant-triad. The role of 
the first oblique mode is similar to that of the plane wave. When the amplitude of the 
first mode becomes 0(a 4r+1 ), the (generalized-)parametric-resonance interaction enhances 
the growth of the second and third modes as in the subharmonic-resonance analysis of this 
paper. The frequency difference between the second and third modes becomes larger as 
the propagation angle of the first mode 6 \ is increased. (This phenomenon may be very 
useful to prove the existence of the critical-layer dynamics in boundary layers and free shear 
flows.) When there exist three pairs of oblique modes of different frequencies in a natural 
flow, a clear A-shaped structure is hard to be observed as many unexcited flow- visualization 
experiments report. 

The (generalized-)parametric-resonance growth of the second and third oblique modes 
allows the phase-locked interaction of Wu & Stewart (1996) to occur when their amplitudes 
become 0(<7 7r / 2+1 ). The phase-locked interaction can take place between the second (or 
third) mode and any other lower- frequency oblique modes of the same phase speed. The 
analysis of Wu & Stewart (1996) indicate that the higher- frequency mode (i.e. the 2nd or 3rd 
mode) can accelerate the growth of the lower-frequency ones. This phase-locked interaction 
occurs earlier than the self-interaction. Both in parametric- resonance and phase-locked 
stages, the first mode grows linearly. When the second and third modes become sufficiently 
large, i.e. <9(o- 3r+1 ), due to the continuous parametric-resonance interaction with the first 
mode, the frequency-detuned self-interaction of this analysis becomes important and all 


oblique modes become singular at a finite downstream position. 

It is shown in Appendix C that higher spanwise-harmonics (A £^ 13 )? higher spanwise- 
periodic mean-flow distortion (AUq 4 } ) and nonlinear mean flow (AUq 4 q) 7 which do not exist 
in the upstream region, are nonlinearly generated by the critical-layer interaction and they 
grow very rapidly to become as large as the base mean flow in the inviscid wall layer at the 
downstream position where the distance from the singularity x s — x is of 0(a r ). The mean 
flow jump across the critical layer was found by Haberman (1972). It is also possible to 
show that the nonlinear interaction in the critical layer generates the critical-layer jumps of 
all higher harmonics and they also become of 0(a) very fast in the inviscid wall layer (cf. 
Smith & Bodonyi 1982). 

The present critical-layer analysis shows that (i) a band of (or bands of a couple of) two- 
dimensional (or small-propagation-angle oblique) modes, which dominate the upstream re- 
gion of a boundary-layer-transition process especially when environmental disturbances are 
relatively small, (ii) bands of oblique modes, whose growth is accelerated by subharmonic- 
parametric-resonance, generalized-parametric-resonance and phase-locked interactions, (iii) 
nonlinearly-generated low-frequency modes, which includes the spanwise-periodic mean-flow 
distortion, (iv) nonlinear mean flow and (v) bands of higher harmonics become very large 
as the instabilities evolve downstream. The growth of these linear and nonlinear modes are 
still governed by the parabolic-type (in the streamwise direction) equations of this analysis. 
The frequency-detuning interaction allows wider bands of modes to interact. 

Near the singularity where x s — x is of 0(a r ), or x s — x is of 0(a“ 1 ), all these instabilities 
become as large as the base mean flow in the inviscid wall layer. The flow in the inviscid 
wall layer will then be fully nonlinear. Goldstein & Lee (1992) and Wu et al (1997) showed 
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that the flow in the next stage will be governed by the unsteady, inviscid, three-dimensional 
triple-deck equations. 


This work was supported by the Acoustics Branch at NASA Lewis Research Center, 
contract number NAS3-98008. 

Appendix A. The 0 ( 1 )- viscosity kernel functions A^, A3, K4 and A' 5 

K 2 = e -X(2G/3+fc)C? C 3 _ 2h s C d( 3 {e- 5 [(Ci+2^/3Kl+(C 1 +C 2 )(2C 2 +C3)C3] Cs ) 

J 0 

+•^0 (Ci 5 C 2 ? C 3 I 0 ? 0 )) — ^c(Ci ? C 2 ? C 3 ) GMCi? C 2 > C 3 ) + C 2 , C 3 IC 1 , 0 ) + C 1 C 2 )}] ? (A 1 ) 

K 3 = e -^( 2 Ci/ 3 +C2)C? [ (Cl + G)Ci 2 + 2 h s C l Cl (Ci + 2 C 2 - 2C 3 ) 

Jo 

-2 h. / Cl dCs {e-X'G (g a (Cu C 2 , Ca) + (Ci, C 2 , Cs|0, C 2 )) - Sc(Ci, C 2 , Ca) (ft(Ci, C 2 , Ca) 

+^a + (Ci,C2,Ca|Ci,0)+ 2 C 2 C 3 - C 1 C 2 )}] , (A 2) 

[ 2 ^ 3 + fc s e~ 2 *< Ci+G/sK? || Cl ( 3 ^ + 2C 2 )(9Ci + 5C 2 ) 

+ lC2^ 6 "(Ci,C2,Ca|0,C2,C3;3) + 2(2Ci +C2)^ 6 + (Ci,C2,Ca|Ci,C2,0;l)}] , (A3) 

Ks = e -(2/3)A[(2C 1+ C2) 3 -2Cf] (2Ci + C 2 ) [Ci (Ci + C 2 ) + h s {3Ci (2Ci + C 2 ) 

+2^ + (Ci,C2,Ca|Ci,0,0;l)}], (A 4) 

where we have put, as in (2.3), 

CiSf- 2 !, C 2 = Xi — x 2 , h s = t sin 2 0, (A 5) 
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with 

^a ± (Cl,C2,C3Ka,Ct) = 4h t d<4e ±X(2Ci/3+Cl+C2 ~ Ca “ C ‘ )C? (Ci + ( 3 ~ ( 4 ) 

Jo 

[1 + 25 ((, - (3) (fe + (3 + (. - a)' 2 ] , (A 6) 

*?(C, , 6, Cla, a, fa *) = r dhe^ : i/3[«.W,K2«c 1 <?T(3/a)(2 C , «,)<>£,] (ft _ & 

Jo 

-? (2Ca + C 2 ) {1 + a (Cl - C 3 ) (3Ci + 2Cfc - Cs + 4Cc) (2Ci - 2Ca + Ci + Cs)}l , (A 7) 

a 


£a = (Cl-C3)(Cl+2C2 + 3C 3 ), (A 8) 

9b = 2 AC 2 C 3 (Ci — C 3 ) (Cl + C 2 + C 3) 2 ? (A 9) 

9c = exp [- A {(Ci + 2C2/3) C 2 2 + (2Ci + 3C 2 + 4C 3 /3) Cf + 2 (Ci + C 2 ) C 2 C 3 }] • (A 10) 

Appendix B. The integrals D p , D s and D m in (5.2) and (5.3) 

D v = dx x ( Xl - l) 2 (x a )- 4 - 2i ^ (2®! - l)- 3+i ^ , (B 1) 


D s = dx 1 J dx 2 (xi - l)(®i + x 2 - 1) 3+ll/ ' < [(xi - l) 2 (xi) 3 1</,J (x 2 ) 3 

+ (x 2 - 1) {xi - 1 - 2r(sin 2 9)(x x - x 2 ) j (xa) -3-1 ^ (x 2 ) _3_1,/ ' J ] , (B 2) 

D m = dx\ jH dx 2 (xj - 1) [2 (xj - l) 2 (xi)- 4 " 2 ^ (x 2 )“ 3 -^ (2x! + x 2 - 2)~ 3+i ^ 

+ (x 2 - 1) (x a + x 2 - 2) (Xi)- 3 -^ (x 2 )- 4 - 2i ^ (xi + 2x 2 - 2)- 3+i ^] . (B 3) 

Appendix C. Critical- layer jumps of higher-order nonlinear modes 

Nonlinear critical-layer equation (1.6.3) indicates that the third-order nonlinear mode 
Ref/^eW 2 ^) appears in the expansion (1.6.15) and the fourth-order ones Ref/^o an< ^ 
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ReU^Je' 42 , among many others, also appear in (1.6.16) (also in equation (6.12) in Lee 
1997a). Here we only consider the single-resonant-triad interaction case (J = \ = 0 and 
T] 0 = 0) and Aj will be replaced by A. 

From the critical-layer analysis, we can show that 


A^i?3 (*) = / Ui%dv = it f dxi f dx 2 A(x 1 )A(x 2 )A*(x 1 + x 2 - x) 

J — oo J — oo J — oo 



e -S<?(2Ci/3+< 2 ) [Ci(2Ci + C 2 ) + (sin 2 0)(. . .) + (sin 4 0)(. . .)] , 

(Cl) 


AUq 4 q(x) = I U$ n dr) = 7T [ dxi\A(xi )\ 2 [l - r(sin 2 0)e _2Jc i /3 | , 

J —oo J —OO L ■* 

(C 2) 


A< } (x) = f i r Jo%dr) = -16tt f dxx f dx 2 f dx 3 f dx 4 
J ~ co J — 00 J—o 0 J —00 J —00 



\^A*(x 2 )A(x z )A(x 4 )A'‘(x 3 + x 4 - Z 2 XC 1 C 3 + •••) + •••]» 

(C3) 

where 



-a 

cow 

II 

i H MZ) c iXo/2 tj( 4 ) H - (4 ) (4) . -„2 f T W 

KT w j\r w cJ ll3 ’ °’° “ 2 K 2 ^^ 0 - 0 ’ Uo,4 ~ *t w cH U 0A , 

(C 4) 


H = 27t Y c a/(kT w cf, 

(C 5) 


c 3 = a : 2 - z 3 , 

(C6) 


Ci and C 2 are defined in (2.3) and M is defined in (1.7.5). 

The nonlinear modes o- 4p+1 Ref7 1 ^e i */ 2 cos3Z, a Sr+1 uf$ and a 5r+1 U^ cos4Z should 
be included in the streamwise velocity expansion (1.10) in order to match with these critical- 
layer jumps. Near the singular point, the amplitude A is given by (5.1). Thus, these higher- 
order velocities become as large as the base mean flow, which is 0(a), in the inviscid wall 
layer when x s — x is 0(a r ). 
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Fig. 1 Quasi-equilibrium interaction between two resonant- triads: In |j4q I > In |^4i j and In |i?o| = 
In |Si| vs. x, curves (i), (ii) and (iii) respectively, when clq = a\ = 0.01 exp(i7r/4), b 0 — h = 1, 
A = 1, k oh = 4/5, k 2 d — 1 and (a) x = 5, (b) x = 10, (c) x = 50 and (d) x = 500 
(dotted, single resonant-triad interaction; dashed, two resonant-triad interaction with x = 0; 
dot-dashed, linear oblique-mode amplitude). 
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Fig. 1 (c) and (d). See previous page for captions. 
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Fig. 2 Quasi-equilibrium interaction between three resonant-triads: ln|Ao|, ln|^i| and \n\A 2 \ 
vs. x , curves (i), (ii) and (iii) respectively, when x = 5, do = di — a? = 0.01 exp(i7r/4), 
b 0 = bi = 62 = 1, A = 1, k 0 b = 4/5 and i<2d — 1 (dotted, single resonant-triad interaction) 
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Fig. 3 Quasi-equilibrium interaction between five resonant-triads: ln|Ao|, ln|Ai| , In 1^421, lnl^l 
and ln[A|| vs. x, curves (i) - (v) respectively, when x — 5 , do = di = <32 = <23 = 64 = 
0.01 exp(i7r/4), 6 0 = 61 = & 2 = 63 = 64 = 1, A = 1, k 0 b = 4/5 and K2d — 1 (dotted, single 
resonant- triad interaction) 
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Fig. 4 Non-equilibrium interaction between two pairs of oblique modes: ln|Ao| and In |Ai| vs. x , 
curves (i) and (ii) respectively, when ao = ai = 1, A = 0 and (a) x = 0, (b) x = 2, (c) x = 5 
and (d) x — 100 (dotted, single-pair interaction; dashed, two-pair interaction with x = 0). 
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Fig. 5 Non-equilibrium interaction between three pairs of oblique modes: In \Aq\ ) In \ A\ \ and \n\A 2 \ 
vs . x, curves (i), (ii) and (iii) respectively, when 5o = 2i = 02 = 1, A = 0 and (a) x = 2 and 
(b) x = 5 (dotted, single-pair interaction; dashed, three-pair interaction with \ — 0). 
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Fig. 7 Non-equilibrium interaction between two pairs of oblique modes: In \Aq\ and In \ Ai\ vs. x } 
curves (i) and (ii) respectively, when x = 2, A = 0, ao = 1 and (a) a\ — 0.8 and (b) a\ — 0.5 
(dotted, single-pair interaction). 
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Fig. 8 Non-equilibrium interaction between two pairs of oblique modes: In |j 4 0 | and In ]Ai\ vs. x, 
curves (i) and (ii) respectively, when \ = 2, A = 0, a a = 1 and (a) a 0 = 0.8 and (b) a 0 = 0.5 
(dotted, single-pair interaction). 
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Fig. 10 Non-equilibrium interaction between three pairs of oblique modes: (a) In |j4i|, (b) In |Ao|, 
(c) \Ao/A\\, (d) In |j 4 2 | and (e) |j4 2 Mi| vs. x when ai = 1, a 0 = a 2 = 0.7, A = 0 and x = 0, 
1, 2, 3 and 5, curves (i) - (v) respectively (dotted, single-pair interaction). 
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Fig. 10 (b) and (c). See previous page for captions. 
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Fig. 10 (d) and (e). See p.70 for captions. 
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Fig. 11 Non-equilibrium interaction between 51 pairs of oblique modes: (a) In \A 3 1 vs. j at x = —2, 
— 1, 0, 1, 1.52, 2, 2.52 and 2.96, curves (i) - (viii) respectively, and (b) ln|j4(_20)|i ln|A(_io)|, 
In |i4o|, In \A\q\ and In I-A 20 I vs. x> curves (i) - (v), when x — 0.2, ni = 0 and A = 0. 
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Fig. 12 Non-equilibrium interaction between 51 pairs of oblique modes: (a) In \A } \ vs. j at x = —2, 
-1, 0, 1, 1.52, 2, 2.52, 3 and 3.32, curves (i) - (ix) respectively, and (b) In \A(_ 2 o)\, In |^4(_io)|, 
In |j4q | i In |j4io| and In |>l 2 o| vs. x, curves (i) - (v), when \ = 0-5, ni = 0 and A = 0. 
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Fig. 13 Non-equilibrium interaction between 51 pairs of oblique modes: (a) In \ A : \ vs. j at x = — 2, 
— 1, 0, 1, 1.52, 2, 2.52 and 2.8, curves (i) - (viii), and (b) In |^4(_io)|, In In \Aq\, In \A$\ 

and In \A\q\ vs. x, curves (i) - (v), when x — 0.2, ni = 5, = —5, C 2 = 1 and A = 0. 
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Fig. 15 Non-equilibrium interaction between 51 pairs of oblique modes: lnjA^ vs. j at x — — 2 
— 1, 0, 1, 1.52, 2, 2.52 and 2.88, curves (i) - (viii) respectively, when x = 0.2, n\ = 0, = 5 

C 2 = 0.8 and A = 0. 
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Fig. 17 Non-equilibrium interaction between three pairs of oblique modes: (a) ln|A/7o 2 jf! an< 
(b) In |A vs. x when a 0 = ai = a 2 = 1, A = 0 and x = 2. ( j,£ ) = (-2,2), (-1,1), (-1,2) 
(0,0), (0,1), (0,2), (1,0), (1,1) and (2,0), curves (i) - (ix) respectively. 
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